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ABSTRACT 

We study two related problems in the context of a supergravity dual to J\f = 1 SYM. One 
of the problems is finding kappa symmetric D5-brane probes in this particular background. 
The other is the use of these probes to add flavors to the gauge theory. We find a rich and 
mathematically appealing structure of the supersymmetric embeddings of a D5-brane probe 
in this background. Besides, we compute the mass spectrum of the low energy excitations 
of M = 1 SQCD (mesons) and match our results with some field theory aspects known from 
the study of supersymmetric gauge theories with a small number of flavors. 
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1 Introduction 



The gauge/string correspondence, an old proposal due to 't Hooft [1], is now well understood 
in the context of maximally supersymmetric super Yang-Mills (SYM) theories. Indeed, the 
so-called AdS/CFT correspondence is a conjectured equivalence between type IIB string 
theory on AdS 5 x S 5 and Af = 4 SYM theory [2]. In the large 't Hooft coupling limit, the 
Af = 4 SYM theory is dual to the type IIB supergravity background corresponding to the 
near-horizon geometry of a stack of parallel D3-branes, whose metric is precisely that of the 
AdSs x S* 5 space. There are nowadays a lot of non-trivial tests of this duality (for a review 
see [3]). 

The extension of the gauge/string correspondence to theories with less supersymmetries 
is obviously of great interest. A possible way to obtain supergravity duals of SYM theories 
with reduced supersymmetry is to consider branes wrapping supersymmetric cycles of Calabi- 
Yau manifolds [4]. In order to preserve some supersymmetry the normal bundle of the cycle 
within the Calabi-Yau space has to be twisted [5]. Gauged supergravities in lower dimensions 
provide the most natural framework to implement this twisting. In these theories the gauge 
field can be used to fiber the cycle in which the brane is wrapped in such a way that some 
supersymmetries are preserved. 

In this paper we will restrict ourselves to the case of the supergravity dual of M = 1 SYM. 
This background, which corresponds to a fivebrane wrapping a two-cycle, was obtained in 
ref. [6] from the solution found in ref. [7] representing non-abelian magnetic monopoles 
in four dimensions. The geometry of this background is smooth and leads to confinement 
and chiral symmetry breaking. Actually, if only the abelian part of the vector field of seven 
dimensional gauged supergravity is excited, one obtains a geometry which is singular at 
the origin and coincides with the smooth one at large distances, i.e. in the UV. Therefore, 
the singularity at the origin is resolved by making the gauge field non-abelian, in complete 
analogy to what happens with the resolution of the Dirac string by the 't Hooft-Polyakov 
monopole. Moreover, as argued in ref. [8], the same mechanism that de- singular izes the 
supergravity solution also gives rise to gaugino condensation. Based on this observation, the 
NSVZ beta function can be reproduced at leading order [9, 10, 11]. Other aspects of this 
supergravity dual have been studied in ref. [12] (for a review see [13]). 

Most of the analysis carried out with the background of [6] do not incorporate quarks in 
the fundamental representation which, in a string theory setup, correspond to open strings. 
In order to introduce an open string sector in a supergravity dual it is quite natural to add D- 
brane probes and see whether one can extract some information about the quark dynamics. 
As usual, if the number of brane probes is much smaller than those of the background, one 
can assume that there is no backreaction of the probe in the bulk geometry. In this paper 
we follow this approach and we will probe with D5-branes the supergravity dual of Af = 1 
SYM. The main technique to determine the supersymmetric brane probe configurations is 
kappa symmetry [14], which tells us that, if e is a Killing spinor of the background, only 
those embeddings for which a certain matrix T K satisfies 

T K e = e. (1.1) 

preserve the supersymmetry of the background [15]. The matrix T K depends on the metric 



1 



induced on the worldvolume of the brane. Therefore, if the Killing spinors e are known, we 
can regard (1.1) as an equation for the embedding of the brane. 

The starting point in our program will be the determination of the Killing spinors of 
the background. It turns out that a simple expression for these spinors can be obtained if 
one considers a frame inspired by the uplifting of the metric from gauged seven dimensional 
supergravity. The realization of the topological twist in this case is similar to the one 
introduced in ref. [16], which generalizes that of ref. [17], to obtain manifolds of G2 holonomy 
and the deformed and resolved conifold from gauged eight-dimensional supergravity. The 
resulting Killing spinors are characterized by a series of projections and all we have to do is 
to find those configurations for which the kappa symmetry condition (1.1) follows from the 
projections satisfied by the Killing spinors. 

The probes we are going to consider are D5-branes wrapped on a two-dimensional sub- 
manifold. We will be able to find some differential equations for the embedding which are, 
in general, quite complicated to solve. The first obvious configuration one should look at is 
that of a fivebrane wrapped at a fixed distance from the origin. In this case the equations 
simplify drastically and we will be able to prove a no-go theorem which states that, unless 
we place the brane at an infinite distance from the origin, the probe breaks supersymmetry. 
This result is consistent with the fact that these M = 1 theories do not have a moduli space. 
In this analysis we will make contact with the two-cycle considered in ref. [10] and show 
that it preserves supersymmetry at an asymptotically large distance from the origin. 

Guided by the negative result obtained when trying to wrap the D5-brane at constant 
distance, we will allow this distance to vary within the two-submanifold of the embedding. To 
simplify the equations that determine the embeddings, we first consider the singular version 
of the background, in which the vector field of the seven dimensional gauged supergravity is 
abelian. This geometry coincides with the non-singular one, in which the vector field is non- 
abelian, at large distances from the origin. By choosing an appropriate set of variables we will 
be able to write the differential equations for the embedding as two pairs of Cauchy-Riemann 
equations which are straightforward to integrate in general. Among all possible solutions 
we will concentrate on some of them characterized by integers, which can be interpreted as 
winding numbers. Generically these solutions have spikes, in which the probe is at infinite 
distance from the origin and, thus, they correspond to fivebranes wrapping a non compact 
submanifold. Moreover, these configurations are worldvolume solitons and we will verify 
that they saturate an energy bound [18]. 

With the insight gained by the analysis of the worldvolume solitons in the abelian back- 
ground we will consider the equations for the embeddings in the non-abelian background. 
In principle any solution for the smooth geometry must coincide in the UV with one of the 
configurations found for the singular metric. This observation will allow us to formulate an 
ansatz to solve the complicated equations arising from kappa symmetry. Actually, in some 
cases, we will be able to find analytical solutions for the embeddings, which behave as those 
found for the singular metric at large distance from the origin and also saturate an energy 
bound, which ensures their stability. 

One of our motivations to study brane probes is to use these results to explore the quark 
sector of the gauge/gravity duality. Actually, it was proposed in refs. [19, 20] that one can 
add flavor to this correspondence by considering spacetime filling branes and looking at their 
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fluctuations. In ref. [21] this program has been made explicit for the AdS§ x S 5 geometry 
of a stack of D3-branes and a D7-brane probe. When the D3-branes of the background and 
the D7-brane of the probe are separated, the fundamental matter arising from the strings 
stretched between them becomes massive and a discrete spectrum of mesons for an Af = 2 
SYM with a matter hypermultiplet can be obtained analytically from the fluctuations of a D7- 
brane probe. In ref. [22] a similar analysis was performed for the Af = 1 Klebanov-Strassler 
background [23], while in refs. [24, 25] the meson spectrum for some non-supersymmetric 
backgrounds was found (for recent related work see refs. [26, 27]). 

It was suggested in ref. [28] that one possible way to add flavor to the Af = 1 SYM 
background is by considering supersymmetric embeddings of D5-branes which wrap a two- 
dimensional submanifold and are spacetime filling. Some of the configurations we will find 
in our kappa symmetry analysis have the right ingredients to be used as flavor branes. They 
are supersymmetric by construction, extend infinitely and have some parameter which deter- 
mines the minimal distance between the brane probe and the origin. This distance should be 
interpreted as the mass scale of the quarks. Moreover, these brane probes capture geomet- 
rically the pattern of R-symmetry breaking of SQCD with few flavors [29]. Consequently, 
we will study the quadratic fluctuations around the static probe configurations found by 
integrating the kappa symmetry equations. We will verify that these fluctuations decay 
exponentially at large distances. However, we will not be able to define a normalizability 
condition which could give rise to a discrete spectrum. The reason for this is the exponential 
blow up of the dilaton at large distances. Actually, this same difficulty was found in ref. 
[30] in the study of the glueball spectrum for this background. As proposed in ref. [30], we 
shall introduce a cut-off and impose boundary conditions which ensure that the fluctuation 
takes place in a region in which the supergravity approximation remains valid. The resulting 
spectrum is discrete and, by using numerical methods, we will be able to determine its form. 

The organization of this paper is the following. In section 2 we introduce the supergravity 
dual of Af = 1 SYM. The Killing spinors for this background are obtained in appendix A, 
where we also obtain those corresponding to the background of refs. [31, 32, 33], which 
represents D5-branes wrapped on a three-cycle. In section 3 we obtain the kappa symmetry 
equations which determine the supersymmetric embeddings. In section 4 we obtain the no-go 
theorem for branes wrapped at fixed distance. In section 5 the kappa symmetry equations 
for the abelian background are integrated in general and some of the particular solutions are 
studied in detail. Section 6 deals with the integration of the equations for the supersymmetric 
embeddings in the full non-abelian background. The spectrum of the quadratic fluctuations 
is analysed in section 7. The asymptotic form of these fluctuations is obtained in appendix 
B. Finally, in section 8 we summarize our results, draw some conclusions and discuss some 
lines of future work. 

1.1 Reader's guide 

Given that this is a long paper, we feel it would be useful to include here a "roadmap" 
to help the reader to quickly find the results of his/her particular interest. Those readers 
interested in the supersymmetry preserved by the background and in the application of kappa 
symmetry to find compact and non compact embeddings in the geometry dual to Af = 1 
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SYM, should pay special attention to sections 2-6 and appendix A. Readers more interested 
in the gravity version of the addition of flavors to J\f = 1 SYM should take for granted section 
3 and look at the solutions exhibited in eqs. (5.19), (5.22) and (6.15), which are what we 
called "abelian and non-abelian unit- winding solutions". Then, they should go straight to 
section 7 and take into account the results of appendix B. 

2 The supergravity dual of AT = 1 Yang-Mills 

The supergravity solution we will be dealing with corresponds to a stack of iV D5-branes 
wrapped on a two-cycle. It can be obtained [6, 7] by considering seven dimensional gauged 
supergravity, which is a consistent truncation of ten dimensional supergravity on a three- 
sphere. To get this background one starts with an ansatz for the seven dimensional metric 
which has a term corresponding to the metric of a two-sphere and looks for a supersymmetric 
solution of the equations of motion. After uplifting to ten dimensions one ends up with a 
solution of type IIB supergravity which preserves four supersymmetries. The ten dimensional 
metric in Einstein frame is: 



ds 2 = e2 



dxl 3 + e 2h ( d6 2 + sin 2 6dip 2 ) + dr 2 + ^ (w l - A 1 ) 2 ] , (2.1) 



where is the dilaton, the unwrapped coordinates x M have been rescaled and all distances 
are measured in units of Ng s a'. The angles G [0, ir] and tp G [0, 2tt) parametrize the 
two-sphere of gauged seven dimensional supergravity. This sphere is fibered in the ten 
dimensional metric by the one-forms A % {% = 1,2,3), which are the components of the non- 
abelian gauge vector field of the seven dimensional supergravity. Their expression can be 
written in terms of a function a(r) and the angles (0, ip) as follows: 

A 1 — -a(r)d0 , A 2 = a(r) sin 6d<p , A 3 = - cos 0dp . (2.2) 

The tu"s appearing in eq. (2.1) are the su(2) left-invariant one-forms, satisfying dw l = 
— | e^k uji Aw k , which parametrize the compactification three-sphere and can be represented 
in terms of three angles <p, and ip: 

w 1 = cos ipd0 + sin ip sin 0d(p , 
w 2 = — sin ipd9 + cos ip sin 0d(p , 

w 3 = dip + cos0dip . (2.3) 

The three angles <p, and ip take values in the rank < (p < 2ir, < < n and < ip < 4n. 
For a metric ansatz such as the one written in (2.1) one obtains a supersymmetric solution 
when the functions a(r), h(r) and the dilaton are: 

2r 

Gin = , : - , 

v ; sinh2r 
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r 2 I 

e 2h = rcoth2r 5 , 

sinh 2r 4 

e - 2 ^ = e -20o^L , (2.4) 
sinn 2r 

where 0o is the value of the dilaton at r = 0. Near the origin r = the function e 2h behaves 
as e 2h ~ r 2 and the metric is non-singular. The solution of the type IIB supergravity includes 
a Ramond-Ramond three-form F( 3 ) given by 

F {3) = -hw l -A l ) A(w 2 -A 2 ) A(w 3 -A 3 ) + - £ F a A(w a -A a ) , (2.5) 

where F a is the field strength of the su(2) gauge field A a , defined as: 

F a = dA a + l -e abc A h AA c . (2.6) 
The different components of F a are: 



F 1 = -a' dr A dO , F 2 = a' sin 9dr A dip , F 3 = ( 1 - a 2 ) sin0d0 A dip , (2.7) 

where the prime denotes derivative with respect to r. Since dF^ = 0, one can represent F( 3 ) 
in terms of a two- form potential C( 2 ) as F( 3 ) = dC^y Actually, it is not difficult to verify 
that C( 2 ) can be taken as: 



C, 



(2) 



■0 ( sin 0C/6 1 A dip — sin ^rf^ A dp) — cos 6* cos Qdip A dp — 



—a ( d6 A w 1 — sin 9dip A w 2 ) 



(2.8) 



Moreover, the equation of motion of F( 3 ) in the Einstein frame is d(^e^*F^ ) = 0, where * 
denotes Hodge duality. Therefore it follows that, at least locally, one must have 

ef*F {3) = dC {6) , (2.9) 

with C(6) being a six-form potential. It is readily checked that C(%) can be taken as: 

C( 6) = dx° A dx l A dx 2 A dx 3 A C , (2.10) 

where C is the following two-form: 

C = \((a 2 -l)a 2 e- 2h - 16e 2 M cosOdpAdr - (a 2 - 1) e~ 2h w 3 A dr + 

g L V / 



+ a' ( sin 6»c?v? Aw 1 + d6 Aw 2 ) 



(2.11) 
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The Killing spinors e of the above background are worked out in Appendix A. They can 
be obtained by requiring the vanishing of the supersymmetry variations of the fermionic 
fields of type IIB supergravity. This requirement leads to a system of first-order BPS dif- 
ferential equations (eq. (A. 16)) for the functions 0, h and a of the ansatz written in eqs. 
(2.1) and (2.5). It can be easily checked that the functions of eq. (2.4) satisfy the system 
(A. 16). Moreover, it follows from the analysis of appendix A that the Killing spinors e are 
characterized by the following set of projections: 

F x () -x 3 ( coso;r 1 2 + sino;r 1 f2)e = e , 

r 12 e = r 12 e , 



e = ie* , (2.12) 

where the T-matrices refer to the frame (A. 2) and the explicit expression of the angle a, 
which depends on the radial coordinate r, is given in eqs. (A. 17) and (A. 25) (see eqs. (A. 3), 
(A. 4) and (A.26)). It is interesting to write here the UV and IR limits of a, namely 

71 

lim a = , lim a = . (2.13) 

r^oo ' r ^0 2 

The BPS equations (A. 16) also admit a solution in which the function a(r) vanishes, i.e. 
in which the one-form A has only one no n- vanishing component, namely A 3 . We will refer 
to this solution as the abelian Af = 1 background. Its explicit form can be easily obtained 
by taking the r — > oo limit of the functions given in eq. (2.4). Notice that, indeed a(r) — > 
as r — > oo in eq. (2.4). Neglecting exponentially suppressed terms, one gets: 

e 2h = r - 1 , (a = 0) , (2.14) 

while <fi can be obtained from the last equation in (2.4). The metric of the abelian background 
is singular at r = 1/4 (the position of the singularity can be moved to r = by a redefinition 
of the radial coordinate). This IR singularity of the abelian background is removed in the 
non-abelian metric by switching on the A 1 , A 2 components of the one-form (2.2). Moreover, 
when a = 0, the angle a appearing in the expression of the Killing spinors (and in the 
projection (2.12)) is zero, as follows from eq. (A. 17). 



3 Kappa symmetry 

As mentioned in the introduction, the kappa symmetry condition for a supersymmetric 
embedding of a D5-brane probe is Y K e = e (see eq. (1.1)), where e is a Killing spinor of the 
background. For e such that e = ie* and when there is no worldvolume gauge field, one has: 

P — — - f m l- m 6 ry (•} Y) 
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where g is the determinant of the induced metric g mn on the worldvolume 

9mn = d m X^d n X u G llv , (3.2) 

with G^y being the ten-dimensional metric and 7 mi ... mg are antisymmetrized products of 
worldvolume Dirac matrices j m , defined as: 

lm = d m X»E v -T v _. (3.3) 

The vierbeins Ej- are the coefficients which relate the one-forms e- of the frame and the 

differentials of the coordinates, i.e. e- = E^dX^ . Let us take as worldvolume coordinates 

" ~ ~ 

(x°, • • • , x 3 , 9, (p). Then, for an embedding with 9 = 9(9, p>), (p = p>(9,(p), tp = ip(9,(p) and 
r = r(9,<p), the kappa symmetry matrix T K takes the form: 

e <t> 

T K = — 1= F x o... x 3 7<^ , (3.4) 

with being the antisymmetrized product of the two induced matrices 70 and 7^, which 
can be written as: 

e-l l9 = e h T x + (V ie + a - ) f ! + V20 f 2 + V 3e f 3 + d e rV r , 

^4j7v> = e h T 2 + + (V 2ip -£)f 2 + ^f 3 + -^r r , (3.5) 

where the Vs can be obtained by computing the pullback on the worldvolume of the left 
invariant one- forms w l , and are given by: 

Vie = ~ cos ip d e 9 + - sin ip sin 9 d g (p , 

1 1 ~ 

sin 9 V iip = - cos ip 8^9 H — sin ip sin 9 d v <p , 

= — \ sin + i cos ^ sin (p , 

sinflV^ = — 7- simp 9 + ^ cos ^ sin 6* d v <p , 

V 3e = ^d e ip + ^cos0d 9 £ , 

sinflV^ = -d^V + ^ cos 9 8^ ip + - cos6> . (3.6) 

By using the projections (A. 3) and (A. 8) one can compute the action of 7^ on the Killing 
spinor e. It is clear that one arrives at an expression of the type: 



e 2 
sin6> 



le^e = [c 12 r 12 + c l2 rif 2 + c^Tifi + c 13 Tif 3 + 

+ c i3 f 13 + c 23 f 23 + c 23 r 2 f 3 ] e , (3.7) 
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where the c's are coefficients that can be explicitly computed. By using eq. (3.7) we can 
obtain the action of T K on e and we can use this result to write the kappa symmetry projection 
T K e = e. Actually eq. (1.1) is automatically satisfied if it reduces to the first equation in eq. 
(2.12). If we want this to happen, all terms except the ones containing r\ 2 e and r\ T 2 e on 
the right-hand side of eq. (3.7) should vanish. Then, we should require 

C li = C 13 = C 13 = C 23 = C 23 = • (3-8) 

By using the explicit expressions of the c's one can obtain from eq. (3.8) five conditions that 
our supersymmetric embeddings must necessarily satisfy. These conditions are: 

e h (V lv + V 2B )=0, (3.9) 

(X d T' 

e h (V 3v , + cos ad 6 r) + (V 2(p - -) smad e r - V 26 sin a -r— = , (3.10) 

2 smd 

(V w + ^)V 3v - V 39 V lip - e h smad e r + 

d d v 

+ ( V 2ip - - ) cosacV - V 2d cosa = , (3.11) 
2 sin 9 

a 

V 3(p V 2S - V 3 e ( V 2(f - - ) - V ltp cos a d e r + 



+ (e h sma + {V w + ^)cosa) ^- =0 , (3.12) 
V 2 J sine* 

e h cosa - (V w + -)sina ) — p — n = . (3.13) 
2 / smO 

Moreover, if we want the kappa symmetry projection T K e = e to coincide with the SUGRA 
projection, the ratio of the coefficients of the terms with Ti T 2 e and r\ 2 e must be tan a, i. e. 
one must have: 



tana = . (3.14) 

Cl2 



The explicit form of C\ 2 and c 12 is: 



G/ T , T r \ a 2 1 t" i — v t cLr 



c i2 = e + VujV^ - V 26 V 1(p - - (V w - V 2(p ) - — - cosaV 3(p d e r -^cosa^^- , 

2 4 smfc> 

d v 

c 12 = e h ( V 2ip - Vw - a ) - sin aV^cV + sm.o>V 3e -r-r . (3.15) 

sin u 

Amazingly, except when r is constant and takes values in the interval < r < 00 (see section 
4), eq. (3.14) is a consequence of eqs. (3.9)-(3.13). Actually, by eliminating V 3 $ of eqs. (3.12) 
and (3.13), and making use of eqs. (3.9) and (3.10), one arrives at the following expression 
of tana: 

e h {V2 V - Vie -a) ,„ 1R , 

tana = T . 3.16 

e 2h + V w V 2ip - V 2e V lip - f (V w - V 2(f ) - % 
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Notice that the terms of (ci 2 ) which do not contain sin a (cos a) are just the ones in the 
numerator (denominator) of the right-hand side of this equation. It follows from this fact 
that eq. (3.14) is satisfied if eqs. (3.9)-(3.13) hold. Moreover, by using the values of cos a 
and sin a given in eq. (A. 17), one obtains the interesting relation: 

(l + a 2 + Ae 2h )(V ie - V 2ip )=Aa[vi e + V ie V 2(p - 1) . (3.17) 

The system of eqs. (3.9)-(3.13) is rather involved and, although it could seem at first sight 
very difficult and even hopeless to solve, we will be able to do it in some particular cases. 
Moreover, it is interesting to notice that, by simple manipulations, one can obtain the fol- 
lowing expressions of the partial derivatives of r: 

d e r = -cosaV^ + sinae~ h [(V w + | ) V 3ip - V 3 oV 1(p ] , 

d v r = cos a sin 9V 3e + sin a sin9e~ h [(V 2ip - ^ ) V 38 - V 3ip V 2e } , (3.18) 
which will be very useful in our analysis. 



4 Branes wrapped at fixed distance 

In this section we will consider the possibility of wrapping the D5-branes at a fixed distance 
r > from the origin. It is clear that, in this case, we have dgr = d^r = and many of 
the terms on the left-hand side of eqs. (3.9)-(3.13) cancel. Moreover e h is non-vanishing 
when r > and it can be factored out in these equations. Thus, the equations (3.9)-(3.13) 
of kappa symmetry when the radial coordinate r is constant and non-zero reduce to: 

V lv + V 2e = V 3<p = V 3e = . (4.1) 

From the equations V 3tp = V 3 e = we obtain the following differential equations for ip 

doip — — cos 9 dgip , dpip — —cos dd^ifi — cos 9 . (4.2) 

The integrability condition for this system gives: 

d v 9d e <p - d e 9d v <p = ^4 • ( 4 - 3 ) 

sm9 

By using this condition and the definition of the V^'s (eq. (3.6)) one can prove that 

V w V 2ip - V 1<P V 26 = ~ . (4.4) 

Let us now define A as follows: 

V 2lf - V w = A . (4.5) 
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By using the expression of the V's in terms of the angles, one can combine eq. (4.5) and the 
condition Vi v + V 2 g = in the following matrix equation 



/ cos ip sin ip 
\ — sin ip cos tp 



( sm.9de9 — smOd^ip 
^ sin 9 sin 9de(p + 8^9 





/-2Asin#\ 


1- 


V o J 



(4.6) 



Since the matrix appearing on the left-hand side is non-singular, we can multiply by its 
inverse. By doing this one arrives at the following equations: 

~ sin 
dad — — — -cL</5 = — 2 A cos^ , 



sin 9 



sin 9 



+ sin9do(p = — 2 A sin^ . 



(4.7) 



Substituting the derivatives of 9 obtained from the above equations into the integrability 
condition (4.3) we obtain after some calculation 



sin 2 9 d e (p + A — 



sin-0 



sin6> 



+ 



sin 2 fl 
sin 2 ^ 



d v (p — A cos^ — 



sin 9 



sin 9 



= A 2 



1 . 



(4.8) 



The right-hand side of eq. (4.8) is non-negative. Then, one obtains a bound for A: 

A 2 > 1 . (4.9) 

Notice that we have not imposed all the requirements of kappa symmetry. Indeed, it still 
remains to check that the ratios between the coefficients c\i and is the one corresponding 
to the projection of the background. Using eq. (4.4) and the definition of A (eq. 4.5)), one 
obtains: 

C12 = e m + a— — , 



C 12 



e h (A - a) 



Then, one must have: 



tana = 



e h (A - a) 



ae 



e 2h + a # 



a 2 +l 



e 2h + 1=0? > 



(4.10) 



(4.11) 



where we have used the values of sin a and cos a given in the appendix A (eq. A. 17). If e h 
is nonzero (and finite), we can factor it out in eq. (4.11) and obtain the following expression 
of A: 

A = — . (4.12) 

1 + a 2 + 4e 2h y ' 

Notice that A depends only on the coordinate r and is a monotonically decreasing function 
such that 0<A<lforO<r<oo and 

1 , 



lim A 

r-»0 



lim A 



. 



(4.13) 
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As A < 1, the bound (4.9) is not satisfied and, thus, there is no solution to our equations 
for < r < oo. Notice that this was to be expected from the lack of moduli space of the 
M = 1 theories. 

Let us now consider the possibility of placing the brane probe at r — > oo. Notice that in 
this case eq. (4.11) is satisfied for any finite value of A. However, the value A = 1 is special 
since, in this case, the right-hand side of eq. (4.8) vanishes and we obtain two equations that 
determine the derivatives of <p, namely: 

a ~ sin ^ a ~ , sin6> 

°e<P = — r^f , d v <p = cosvp (4.14) 

sin 9 smO 

Using these equations into the system (4.7) for A = 1 one gets the following equations for 
the derivatives of 9: 

dg9 = — cos-0 , dj) = — sin ^ sin , (4-15) 

and, similarly, the equations (4.2) for ip become: 

dgip = sin^cot^ , dtpt/j = — sin 9 cot 9 cos ip — cos9 (4.16) 

The equations (4.14) and (4.15) can be regarded as coming from the following identifications 
of the frame forms in the (9, ip) and (9, <p>) spheres: 



( d6 \ ( cos ip — sin ip \ ( —d9 



\ sin 9d(p ) \ sin ip cos ip ) \ sin 6 dip t 



(4.17) 



The differential equations (4.16) are just the integrability conditions of the system (4.17). 
Another interesting observation is that one can prove by using the differential eqs. (4.14)- 
(4.16) that the pullbacks of the su(2) left-invariant one- forms are 

P^ 1 ] = -d6 , P[w 2 } = sin9dip , P[w 3 } = - cos 9dp . (4.18) 

Let us try to find a solution of the differential equations (4.14)-(4.16) in which 9 = 9(9) 
and <p = <p(ip). The vanishing of d v 9 and de<p immediately leads to sin^ = or ip = 
0,7r (mod 2ir). Thus ip is constant in this case. Let us put cos^ — r] — ±1. The vanishing 
of deip is automatic, whereas the condition d v ip = leads to a relation between 9 and 9: 

cot 9 = -T]cot9 (4.19) 

In the case ip — 0, one has rj — 1 and the previous relation yields 9 = n — 9. Notice that this 
relation is in agreement with the first equation in eq. (4.15). Moreover, the second equation 
in (4.14) gives ip — ip. Similarly one can solve the equations for ip — ir. The solutions in 
these two cases are just the ones used in ref. [10] in the calculation of the beta function 
(with some correction in the ip — case to have the correct range of 9 and 9), namely: 

9 = 7r — 9 , <p — ip , ip — (mod 2ir) , 

9 = 9, ip = 2n-p , ^ = vr (mod 2vr) (4.20) 
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It follows from our results that the embedding of ref. [10] is only supersymmetric asymp- 
totically when r — > oo. In this sense, although it is somehow distinguished, it is not unique 
since for any embedding such that the Vs are finite when r — > oo, the determinant of the 
induced metric diverges as y/—g ~ e~ +2h and the only term which survives in the equation 
T K e = e is the one with the matrix Ti 2 , giving rise to the same projection as the background 
for r — > oo. 



5 Worldvolume solitons (abelian case) 

Let us consider the case a = a = in the general equations of section 3. From equations 
(3.9) and (3.17) we get the following (Cauchy-Riemann like) equations: 

V ie = V 2ip , V llf) = -V 2e , (5.1) 

whereas, from eq. (3.18) we obtain that the derivatives of r are given by: 

re = -V 3ip , r v = sin 9V 36 , (5.2) 

where rg = dgr and r v = d 9 r. It can be easily demonstrated that, in this abelian case, 
the full set of equations (3.9)-(3.13) collapses to the two pairs of equations (5.1) and (5.2). 
Notice that c 12 = when a = a = and eq. (5.1) holds and, thus, eq. (3.14) is satisfied 
identically. 

Let us study first the two equations (5.1). By using the same technique as the one 
employed in section 4 to derive eq. (4.7), it can be shown easily that they can be written as 

sm9d e 9 = sm.9d v (p , dj) = —smOsmOdgip . (5.3) 

In order to find the general solution of eq. (5.3), let us introduce a new set of variables u 
and u as follows: 

9 9 

u = log tan - , u = log tan - . (5.4) 

Then, eq. (5.3) can be written as the Cauchy-Riemann equations in the (u, (p) and (u, <p) 
variables, namely: 

du d(p du d(p 

du ' dtp du 

Since u,u G (— oo, +oo) and <p, <p G (0, 2ir), the above equations are the Cauchy-Riemann 
equations in a band. The general solution of these equations is of the form: 

u + tip = f(u + iip) , (5.6) 

where / is an arbitrary function. Given any function /, it is clear that the above equation 
provides the general solution 9(9, (p) and p(9,p) of the system (5.3). 
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Let us turn now to the analysis of the system of equations (5.2), which determines the 
radial coordinate r. By using the explicit values of V 3(p and V30, these equations can be 
written as: 



1 a 1 I cos6 1 

'« = -2^e d ^-2l^e d ^-r oid > 



sin# sin# ~ _ 
r> = ^-^e^H — cos 6» (5.7) 

where 9(9, tp) and tp(9, tp) are solutions of eq. (5.3). In terms of the derivatives with respect 
to variable u defined above (sin#<9,9 = d u ), these equations become: 

r u = -- d^tp - i cos 9 d v <p - - cos6» , 



r v = i d u tp + i cos 9 d u tp . (5.8) 

The integrability condition of these equations is just d v r u = d u r v . As any solution (9, <p) of 
the Cauchy-Riemann equations (5.3) satisfies: 

djidytp = -d u 9d u (p , (5.9) 

and, since (p, being a solution of the Cauchy-Riemann equations, is harmonic in (u, tp), it 
follows that d v r u = d u r v if and only if tp is also harmonic in (u,tp), i.e. the differential 
equation for tp is just the Laplace equation in the (u, tp) plane, namely 

dltp + d 2 u tp = . (5.10) 

Remarkably, the form of r(9, tp) can be obtained in general. Let us define: 

k(9,tp) = / dtpsm9d e tp(9,tp) - / — 0^(9,0) , (5.11) 
j J sin w 

It follows from this definition and the fact that tp is harmonic in (u, tp) that tp and A also 
satisfy the Cauchy-Riemann equations: 

OA = dtp dA = _dtp 

dtp du ' du dtp 

Thus tp and A are conjugate harmonic functions, i.e. tp + iA is an analytic function of u + itp. 
Notice that given A one can obtain tp by integrating the previous differential equations. It 
can be checked by using the Cauchy-Riemann equations that the derivatives of r, as given 
by the right hand side of eq. (5.8), can be written as re = dgF, r v = d v F, where: 



F(9,p) = \ 



A(9,tp) - log (sin 9 sin 9(9, tp)) 



(5.13) 



Therefore, it follows that: 

e 2r = C ; (5. 14 ) 

sin 9 sin 9(9, tp) 

with C being a constant. We will make use of this amazingly simple expression to derive the 
equation of some particularly interesting embeddings. 



13 



-6 -4 -2 2 4 6 x 

Figure 1: Curves y = y(x) for three values of the winding number n: n = (solid line), 
n = 1 (dashed curve) and n = 2 (dotted line). These three curves correspond to r* = 1. 

5.1 n- Winding solitons 

First of all, let us consider the particular class of solutions of the Cauchy-Riemann eqs. (5.5): 

u + i<p = n{u + i<p) + constant , (5.15) 

where n is an integer and the constant is complex. In terms of the original variables: 

9 ( 9 \ n 
tan- = c ( tan- J , (p — nip + ip , (5.16) 

with c and <po being constants. It is clear that in this solution the (p coordinate of the probe 
wraps n times the [0, 27r] interval as ip varies between and 2n. Let us now assume that the 
coordinate ip is constant, i.e. ip = ipo. It is clear from its definition that the function A(#, ip) 
is zero in this case. Moreover, by using the identities 

2tanf x ll — cos a: . 

siM = ^— , tan— = \ , 5.17 

1 + tan 2 f ' 2 V 1 + cosa; ' v ; 

one can prove that: 

S m9 = 2V~c- J Sm6) " (5.18) 

(1 + cos^)™ + c(l - cos6») n v ' 

where c = c 2 . After plugging this result in eq. (5.14), one obtains the explicit form of the 
function r(9), namely: 

2r e 2r * (1 + cosfl)" + c(l - cos^)" 
6 ~ 1 + c (shitf)^ 1 ' { ' 

where r* = r{9 = ir/2). We will call n-winding embedding to the brane configuration 
corresponding to eqs. (5.16) and (5.19) for a constant value of the angle ip. 

Let us pause for a moment to study the function (5.19). First of all it is easy to verify 
that this function is invariant if we change n — > — n and c — > 1/c (or equivalently changing 
9 — > 7r — 9 for the same constant c). Actually, in what follows we shall take the integration 
constant c = 1 and thus we can restrict ourselves to the case in which n is non-negative. 
In this c = 1 case r* is the minimal separation between the brane probe and the origin. 
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Another observation is that r diverges for 9 = 0,7r, which corresponds to the location of 
the spikes of the worldvolume solitons. Therefore the supersymmetric embedding we have 
found is non-compact. Actually, it has the topology of a cylinder whose compact direction 
is parametrized by <f. This cylinder connects the two poles at 9 = 0, n of the (9, if) sphere 
at r = oo and passes at a distance r* from the origin. 

It is also interesting to discuss the symmetries of our solutions. Recall that the angle ip is 
constant for our embeddings. Thus, it is clear that one can shift it by an arbitrary constant 
e as ip — > ip + e. This U(l) symmetry corresponds to an isometry of the abelian background 
which quantum-mechanically is broken to TL 2 n a s a consequence of the flux quantization 
of the RR two-form potential [6, 34, 35]. In the gauge theory side this isometry has been 
identified [6, 34, 35] with the U(l) R-symmetry of the M = 1 SYM theory, which is broken 
down to TL 2N by a field theory anomaly [29] . On the other hand, it is also clear that we have 
an additional U(l) associated to constant shifts in <f, which are equivalent to a redefinition 
of ipo in eq. (5.16). 

To visualize the shape of the brane in these solutions it is rather convenient to introduce 
the following cartesian coordinates x and y 

x = rcosO , y = rsin9 . (5.20) 

In terms of (x,y) the D5-brane embedding will be described by means of a curve y = y(x). 
Notice that y > 0, whereas — oo < x < +oo. The value of the coordinate y at x = is just 
r*, i.e. y(x = 0) = r*. Moreover, for large values of the coordinate x, the function y(x) — > 
exponentially as 

y(x) w C \x\ e^wr' 1 ' , (\x\ -> oo) , (5.21) 

where C is a constant. To illustrate this behaviour we have plotted in figure 1 the curves 
y(x) for three different values of the winding n and the same value of r*. 

A particularly interesting case is obtained when n — ±1. By adjusting properly the 
constant (po in eq. (5.16) the angular embedding reduces to: 

9 = 9, <p = (f + constant , (n — 1) , 

9 = it — 9 , (p = 2n — if + constant , (n = —1) , (5.22) 

with ip being constant. These types of angular embeddings are similar to the ones considered 
in ref. [10] (although they are not the same, see eq. (4.20) ) and we will refer to them as 
unit-winding embeddings. Notice that the two cases displayed in eq. (5.22) represent the 
two possible identifications of the (9, if) and (9, <f) two-spheres. 

When n = the brane is wrapping the (9, if) sphere at constant values of 9 and <f, i.e. 
one has: 

9 = constant = 9q , if — constant = fo , (n = 0) (5.23) 

We will refer to this case as zero-winding embedding [9]. 

One can verify that the brane embeddings we have found are solutions of the probe 
equations of motion. Actually, they are supersymmetric worldvolume solitons of the D5 
brane probe. To illustrate this fact let us show that these configurations saturate a BPS 



15 



energy bound. To simplify matters, let us assume that the angular embedding is the one 
displayed in eq. (5.16) and let r(9) be an arbitrary function. The Dirac-Born-Infeld (DBI) 
lagrangian density for a D5-brane with unit tension is: 



C = 



a 



(6) 



(5.24) 



where g st is the determinant of the induced metric in the string frame (g st = e 3 ^ g) and 
P C(6) is the pullback on the worldvolume of the RR six-form of the background. The 
elements of the induced metric for the n-winding solution along the angular coordinates are: 

n 2 sin 2 9 



9ee 



6 2 [ e 2h + — 



+ n 



4 sin 2 
n 2 sin 2 6 



(5.25) 



From this expression one immediately obtains the determinant of the induced metric, namely: 



= e 2 smt) 



\ 



o2h 



+ 



n 2 sin 2 9 



+ vi 



4 sin 2 ^ ' Sip J V" ' 4 sin 2 # 



D 2h 



+ 



n 2 sin 2 9 



Moreover, the pullback on the worldvolume of the two- form C is 1 : 



C 



lQe 2h cos 9 - ne~ 2h cos 9)r e d(pAd9 



The hamiltonian density H for a static configuration is just 7i = — C or: 
H = e 2 * 



sin 9 



1 



e 2h + _ 



n 2 sin 2 9 2 
4 sin p 



n 2 sin 2 ^ 



+ — 



4 sin 2 ^ 



(5.26) 



(5.27) 



-(16e 2h cos 9 - ne- 2h cos9)r e 
8 



It can be checked that, for an arbitrary function r(9), one can write TC as: 

H = Z + S , 



(5.28) 



(5.29) 



1 It is worth mentioning that the pullback of the RR two-form to the worldvolume is 

P[C (2) ] = ^dtpA (nsmOS - sin9d9) , 

where 9(6) is the function displayed in cq. (5.18) . From this expression it is straightforward to verify that 
the RR two-form flux through the two-submanifold where we are wrapping our brane is 



/ 



P[C {2) ] = thMM-I) 



and thus it vanishes iff n = ±1. 



16 



where Z is a total derivative: 



-d e 



& [e 2h cos 6 + j cos 



and S is non-negative: 



S > 



(5.30) 



(5.31) 



with S = precisely when the BPS equations for the embedding are satisfied. The expression 
of S is: 



sin Be 2 * 



\ 



2h 



n 2 sin 2 



+ — 



4 sin 2 # 



2h 



+ 



n 2 sin 2 6 1 
T sin 2 



- e 



„2/i 



+ — 



2 sin 2 # 



4 sin 2 # 



- V* 



3tp ' 



(5.32) 



The BPS equation for r in this case is rg = —V^ (see eq. (5.2)). If this equation is satisfied, 
the first term on the right-hand side of eq. (5.32) is a square root of a perfect square which 
cancels against the second term of this equation. Moreover, it is easy to check that the 
condition S > is equivalent to: 



(r + V^f > , 



(5.33) 



which is obviously satisfied and reduces to an equality if and only if the BPS equation for 
the embedding is satisfied. 



5.2 (n,m)-Winding solitons 

The solutions found in the previous section are easily generalized if we allow the angle ip to 
wind a certain number of times as the coordinate tp varies from (p = to ip = 2tt. Recalling 
that ip ranges from to 4n, let us write the following ansatz for ip(p): 



ip — ip + 2mip , 



(5.34) 



where m is an integer. It is obvious that the above function satisfies the Laplace equation 
(5.10). Moreover, its harmonic conjugate A is immediately obtained by solving the Cauchy- 
Riemann differential equations (5.12), namely: 



A = —2mu . 

In terms of the angle 9, the above equation becomes: 

1 



(tanf) 



2/ii 



(5.35) 



(5.36) 
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By plugging this result in eq. (5.14), and using the value of sin 9 given in eq. (5.18), it is 
straightforward to obtain the function r{9) of the embedding. One gets: 

2r (1 + COsfl)" + C(l ~ COS0)" 

6 " 1 + c [tan|] 2m (sin^)«+i ' 1 J 

where, as in the n-winding case, r* = r{6 = n/2). 

An interesting observation concerning the solution we have just found is that, by choosing 
appropriately the winding number to, one of the spikes of the m = solutions at 6 = or 
6 = 7r disappears. Indeed if, for example, n is nonnegative and we take 2m — n + 1, the 
function r{6) is regular at 6 = 0. Similarly, also when n > 0, one can eliminate the spike at 
6 = 7r by choosing 2m = — n — 1. 

5.3 Spiral solitons 

By considering more general solutions of the Cauchy-Riemann equations (5.5) and (5.12) we 
can obtain many more classes of supersymmetric configurations of the brane probe. One of 
the questions one can address is whether or not one can have embeddings in which r is finite 
for all values of the angles. We will now see that the answer to this question is yes, although 
the corresponding embeddings seem not to be very interesting. To illustrate this point, let 
us see how we can find functions ip and A such that they make the radial coordinate of the 
n-winding embedding finite at 6 = 0, 7r. First of all, notice that, in terms of the Cauchy- 
Riemann variables u and u defined in eq. (5.4), we have to explore the behaviour of the 
embedding at u, u — > ±oo. Since 

2e u ~ 2e" 
sm6 = , sinfl = , (5.38) 

one has that sin# — > e - '"', sin^ — > e - '"' as u,u — > ±oo. Then, the factors multiplying e A 
in eq. (5.14) diverge as e'"' + '"L In the n-winding solution \u\ = \n\\u\ and, therefore this 
divergence is of the type e^ n l +1 ^ u L We can cancel this divergence by adding a A such that 
e A — > as u — > ±oo in such a way that, for example, A + (|n| + l)\u\ — > — oo. This is clearly 
achieved by taking a function A such that A — > — u 2 . It is straightforward to find an analytic 
function in the (u, (p) plane such that its imaginary part behaves as — u 2 for u — > ±oo. One 
can take 

iP + iA = -i(u + tip) 2 = 2uip - i(u 2 - ip 2 ) . (5.39) 
From this equation we can read the functions ip and A. In terms of and ip they are: 

6 6 
i[j = 2up = 2p logtan - , A = -u 2 + p 2 = -( logtan-) 2 + p 2 . (5.40) 

In this case r — > 0, ip — > ±oo as 6 — > 0, 7r, which means that we describe an infinite spiral 
which winds infinitely in the ip direction. Notice that, although r is always finite, the volume 
of the two-submanifold is infinite due to this infinite winding. One can try other alternatives 
to make the radial coordinate finite. In all the ones we have analyzed one obtains the infinite 
spiral behaviour described above. 
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6 Worldvolume solitons (nonabelian case) 



Let us consider the full nonabelian background and let us try to obtain solutions to the 
kappa symmetry equations (3.9)-(3.13). Actually we will restrict ourselves to the situations 
in which r only depends on the angle 9. It can be easily checked that, in this case, only four 
of the five equations (3.9)-(3.13) are independent. As an independent set of equations we 
will choose eqs. (3.9), (3.17) and 



d e r 



e h V ? 



3<p 



e h cos a + ( V%p — § ) sin a ' 
sin aV Xip d 6 r - e h V 3e = , 



(6.1) 
(6.2) 



which can be obtained from eqs. (3.10) and (3.13) after taking d^r = 0. 

We will now try to find the non-abelian version of the solutions found in the abelian 
theory for arbitrary winding n. With this purpose, let us consider the following ansatz for 

<p(9,<p) =n V + f(9) , (6.3) 

while we shall assume that 9, ip and r are functions of 9 only. We will require that, in the 
asymptotic UV, <p — > nip. It this clear that in this ansatz d^Cp = n and that dgip = dgf. 
Moreover, from eq. (3.9) we can obtain the relation between dgip and dg9, namely: 



de<p = tanip 



8,9 



n 



sin 9 sin 9 

Using this value of detp, we get the following values of the V functions: 



(6.4) 



V w = 



d ft 9 



Vi 



cos ip 
n sin 9 



— n ■ 



sin 9 sin 2 tp 



2 sin# 
n sin 9 



sin 9 cos ip 
sin^ = -V 2 e , 



2y? 



2 sin# 



cos^ , 



V, 



3(9 



1 1 ~ ~ TX COS 

- dgib + - cot 9 tan ib dg9 — — tan ib - 

2 2 2 sin 9 



= -r 



n cos 9 1 

+ - cot 9 . 

2 sm9 2 



(6.5) 



By using these values in eq. (3.17), one gets the value of dg9 in terms of the other variables: 

n sin 9 cosh 2r — sin 9 cos ib 



d ft 9 



sin 9 cosh 2r — n sin 9 cos ip 



(6.6) 
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On the other hand, by combining the two equations (6.4) and (6.6), we obtain: 



n 2 sin 2 9 - sin 2 9 

Og(f> = ~ smip . (6.7) 

sin 9 sin 9 ( sin 9 cosh 2r — n sin 9 cos ip ) 

Moreover, plugging the values of Vi^, V 2ip , V 3 g and V 3ip in eqs. (6.1) and (6.2) we can obtain 
the values of the derivatives of ip and r. The result is: 

_ . n cot 9 sin 9 + cot 9 sin # . , 

cfoV = sin?/; , 

sin # cosh 2r — n sin # cos tp 

1 n cos # + cos 6* . , „ ,„ 
<9 e r = -- = sinh2r. (6.8) 

2 sin cosh 2r — n sin # cos ^ 

It follows from eq. (6.6) that, asymptotically in the UV, sin 9 dg9 — > nsin#. In order 
to fulfill this relation for arbitrary r it is easy to see from eq. (6.6) that one must have 
(nsin(9) 2 = sin 2 9, which only happens for n = ±1 and sin# = sin#. Noticing that for these 
values one has d$9 = ±1, one is finally led to the two possibilities of eq. (5.22): 9 = 9 for 
n = 1 and 9 = n — 9 for n = —1. Notice that in the two cases of eq. (5.22) this equation 
implies that dg (p = and thus when n — ±1 the angular identifications of the abelian unit- 
winding embeddings (eq. (5.22)) also solve the non-abelian equations (6.6) and (6.7) for all 
r. 

For a general value of n one has that asymptotically in the UV dgtp — > and dgip — > 0, 
as in the abelian solutions. Moreover it follows from eqs. (6.7) and (6.8) that (p and ip can 
be kept constant for all r if sin^ = 0, i.e. when = 0,ir mod 2tt. For this values of ip the 
equations simplify and, although we will not attempt to do it here, one could try to integrate 
numerically the equations of 9 and r. It is however interesting to point out that, contrary to 
what happens in the abelian n-winding solution, the angle ip cannot be an arbitrary constant 
for the non-abelian probes. As we will argue below, this is a geometrical realization of the 
breaking of the R-symmetry of the corresponding M = 1 SYM theory in the IR. On the 
contrary, the angle Cp can take an arbitrary constant value, as in the abelian solution. 



6.1 Non abelian unit-winding solutions 

Let us now obtain the non-abelian generalization of the unit-winding solutions. First of all 
we define 

rj = n = ±1 . (6.9) 

We have already noticed that for unit-winding embeddings the values of 9 and <p displayed in 
eq. (5.22) solve the non-abelian differential equations (6.6) and (6.7). Therefore, let us try 
to find a solution in the nonabelian theory in which the embedding of the (0, <p) coordinates 
is the same as in the abelian theory, i.e. as in eq. (5.22). For this type of embeddings 
sin^ = sin#, d e 9 = rj and eq. (6.5) reduces to: 

rjcosip r/smip 
Vie = — o — ' Vl ? = — o — ' 
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V? e + K = - , V 2 e + V 2 % = - , (6.11) 



f] sin ip rjcostp 

V26 = , V 2(p = , 

^ = y, V 3v = cot6, (6.10) 

where we have denoted ipg = dgip. As a check, notice that Vi$, Vi v , V20 and satisfy eqs. 
(5.1). It follows from eq. (3.17) that they must also satisfy 

I v 2 4- V 2 — — 

which indeed they verify. Moreover, by substituting sin# = sin 6* and cos 6* = rj cos 9 in eq. 
(6.8), we obtain the following differential equations for ip(9) and r(9): 

2i] sin tp 
sinh 2r ' 

cot^ . , . . 

r e = 7 smh2r • 6 - 12 

coshzr — rjcosip 
These equations can be integrated with the result: 




A cothr , 
sinh r C 



V / A 2 + tanh 2 r sin ^ 



(6.13) 



where A and C are constants of integration. Eq. (6.13), together with eq. (5.22), determines 
the unit-winding embeddings of the probe in the non-abelian background. Notice that, as 
in the corresponding abelian solution, r diverges when 9 = 0, ir, i.e. the brane probe extends 
infinitely in the radial direction. On the other hand, it is also instructive to explore the 
r — > 00 limit of the solution (6.13). First of all, it is clear that when r — > 00 the angle i/j 
reaches asymptotically a constant value ipo, given by 

I - A 2 

cos ^0 = 1 + A2 V ■ (6-14) 

Moreover, when r — > 00 the function r{9) displayed in eq. (6.13) becomes, after a proper 
identification of the integration constants, exactly the one written in eq. (5.19) for n — ±1 
and c = 1. Notice that the angle ip in the embedding (6.13) is not constant in general. 
Actually, only when A = or A = 00 the coordinate ip remains constant and equal to 
0,7r mod27r (cosijj = f] for A = and cosV' = — V in the case A = 00). It is interesting to 
write the dependence of r on 9 in these two particular cases. When A = 00 the solution is 



(9, if 77 = +1 , f 
U-9, if 77 = -1 , ' ^ l: 



9, if 77 = +1 , / V 9 + constant, if 77 = +1 , 

2n — ip + constant, if 77 = — 1 



. 7r, 37r, if 77 = +1 , . , sinhr* 

^ = \ n o 1 ' smhr = • a ' ( 6 - 15 ) 

1 0, 27r, if 77 = -1 , sin0 
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Figure 2: Comparison between the non-abelian (solid line) and abelian (dashed line) unit- 
winding embeddings for the same value of r* . The non-abelian embedding is the one 
corresponding to eq. (6.15) and the abelian one is that given in eq. (5.19) for n — 1 and 
c = 1. The curves for two different values of r* (r* = 0.5 and r* = 1) are shown. The 
variables (x,y) are the ones defined in eq. (5.20). 

where r* is the minimal value of r (i.e. r* = r{9 = 7r/2)) and we have also displayed the 
angular part of the embedding. Notice that, for a given sign of the winding number 77, only 
two values of ip are possible. Thus, in this solution, the U(l) symmetry of shifts in ip is 
broken to a S 2 symmetry. This will be interpreted in section 7 as the realization, at the 
level of the brane probe, of the R-symmetry breaking of the gauge theory. 

To have a better understanding of the solution (6.15) we have plotted it in figure 2 in 
terms of the variables (x,y) defined in eq. (5.20). For comparison we have also plotted the 
abelian solution corresponding to the same value of r*. In this figure the embeddings for 
two different values of the minimal radial distance r* are shown. When r* is large enough 
(f* > 2) the two curves become practically identical. 

Let us now have a look at the case of the A = embeddings. The function r(0) in this 
case can be read from eq. (6.13), namely: 

C 

coshr = . (6.16) 

sin 9 

We have plotted in figure 3 the profile for these embeddings in terms of the variables (x, y) 
of eq. (5.20). Notice that, when C is in the interval (l,oo) it can be parametrized as 
C = coshr*, with r* > being the minimal radial distance between the probe and the 
origin. On the contrary, when C lies in the interval [0, 1] the brane reaches the origin when 
sin 9 = C. We have thus, in this case, a one-parameter family of configurations which pass 
through the origin. 

As in their abelian counterparts, these worldvolume solitons for the non-abelian back- 
ground saturate an energy bound. In order to prove this fact, let us define 

D = coth2r - r] . (6.17) 
sinh2r v 1 

Notice that D > for any real ip and r. Moreover the equation for r(9) can be written as 
r e = —cot9/D. For arbitrary functions r(9) and ip(9) the hamiltonian density takes the 
form: 



H = e 2<t> sm9 



7 \2 / cot#\2 rD / 2r]siml) y 

(r-recatO) +rD(r e + —) + ~ O - r) sinh2r cot9 ) 
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Figure 3: Graphic representation of the unit winding embedding of eq. (6.16) for three values 
of the constant C: C = 0.5 (dashed line), C — 1 (solid line) and C = 1.5 (dotted line). The 
variables (x,y) are the ones defined in eq. (5.20). 



-(2e 2h - 1 {a 2 -lfe- 2h ) cot 9 r e 



8 

It can be verified that H can be written as H = Z + S, where 

d 



Z = — 



de 



e 2 * r cos 9 



(6.18) 



(6.19) 



(this is the same value as in the abelian soliton for n — 1). The expression of S is: 
S = e 2(f, sm0 



I / \2 / cot# \2 rD 

(r - r e cot 9) + rD (r e + — — J + — 



2r/ sin ip 



D sinh 2r 



— ^ r — cot ^ ) 



cot ) - 



(6.20) 



As in the abelian case, if the BPS equations (6.12) are satisfied, the square root on eq. (6.20) 
can be exactly evaluated and S vanishes. Furthermore, one can easily check that S > is 
equivalent to the condition 



/ cot 9 \ 2 rD / 
rD [re + ^r- ) ~ —\^e 



D 



4 



2ri sin ib \2 

^' - r cote) > o , 

D sinh 2r ' 



(6.21) 



which, since -D > 0, is trivially satisfied for any functions r(9) and ip(9). Moreover, as 
r — r e cot9 > for the solution of the BPS equations, it follows that our BPS embedding 
saturates the bound. 



6.2 Non abelian zero-winding solutions 

The differential equations for the non-abelian version of the zero-winding solution can be 
obtained by putting n = in our general equations. Actually, by taking n = in the second 
equation in (6.8) one obtains the differential equation which determines the dependence of 
r on the angle 9, namely: 

cot 9 , , 

T > = "2^ • < 6 - 22 > 
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which can be easily integrated, namely: 



sinh 2r 



C 
sin 6* 



(6.23) 



Notice that, as in the abelian case, this solution has two spikes at 9 — 0, tc, where r diverges 
and, thus, the brane probe also extends infinitely in the radial direction. Moreover, the 
minimal value of the radial coordinate, which we will denote by r*, is reached at 9 = tt/2. 
This minimal value is related to the constant C in eq. (6.23), namely sinh 2?% = C. It is 
readily verified that for large r this solution behaves exactly as the zero-winding solution in 
the abelian theory. Moreover, it follows from eq. (6.4) that, in this n = case, the angle p 
only depends on 9. Actually, the differential equations for the angles 9, <p and ip as functions 
of 9 are easily obtained from eqs. (6.6), (6.7) and (6.4): 



d ff 9 



d e p = 



d 6 ip 



cos ip 
cosh2r 

1 sin^ 
cosh2r sin 9 

cot 9 sin ip 
cosh 2r 



(6.24) 



By combining the equations of ip and 9 one can easily get the relation between these two 
angles, namely: 



sin-?/' 



B 
sin 9 



(6.25) 



Notice that, for consistency, B < 1 and sin# > B. We can also obtain p = p(9) and 
9 = 9(9): 



p = — arctan 



Bcos9 



sin 2 9 - B 2 



+ constant , 



arcsm 



cos^ 
VI -B 2 



arcsm 



cos 9 



+ constant 



(6.26) 



Actually, much simpler equations for the embedding are obtained if one considers the par- 
ticular case in which the angle ip is constant. Notice that, as was pointed out after eq. (6.8), 
this only can happen if — 0, 7r (mod 2n) (see also the last equation in (6.24)). These solu- 
tions correspond to taking the constant B equal to zero in eq. (6.25). Moreover, it follows 
from the eq. (6.24) that p is an arbitrary constant in this case, while the dependence of 9 on 
9 can be obtained by combining eqs. (6.23) and (6.24). If we denote cosip = e, with e = ±1, 
one has: 

sinh 2?% . ,~ ~ s cos# ,„ 

(6.27) 



sinh 2r 



sin^ 



sin( 9 — 9* 



cosh 2r* 



where 9* = 9(9 = tt/2). Notice that there are four possible values of ip in this zero-winding 
solution and, thus, the U(l) R-symmetry is broken to 2Z 4 in this case. 
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Let us finally point out that, also in this case, the hamiltonian density H can be put as 
H = Z + S, where S is nonnegative (S = for the BPS solution) and Z is given by: 



e 2<t> cos 9 [ r coth 2r H ^ — 

4 2 sinh 2 2r 



(6.28) 



6.3 Cylinder solutions 

We shall now show that there exists a general class of supersymmetric embeddings for the 
non-abelian background. For convenience, let us consider r as worldvolume coordinate and 
let us assume that the D5-brane is sitting at the north poles of the (6, ip) and (6, (p) two- 
spheres, i.e. at 6 — 6 — 0. In the remaining angular coordinates ip, <p and ip, the embedding 
is characterized by the equation: 

p q s 

where (</?o, <fio, V'o) an d (p,q,s) are constants. Notice that, if one of the constants of the 
denominator in (6.29) is zero, then the corresponding angle must be a constant. Let us 
parametrize these embeddings by means of two worldvolume coordinates a\ and a 2 , defined 
as follows: 

<7i = = = , a 2 — r . (6.30) 

p q s 

It is straightforward to demonstrate that the pullback to the worldvolume of the forms w l 
and A 1 is given by: 

Plw 1 } = P[w 2 } = , P[w 3 } = (q + s)da 1 , 

PIA 1 ] = P[A 2 ] = , P[A 3 ] = -pdo x . (6.31) 

It follows from these results that the pullback of the frame one-forms e l and e J is zero for 
i,j = 1,2, whereas P[e 3 ] is non-vanishing. As a consequence, the induced Dirac matrices 
are: 

l.i = ^(P + q + s) ^ f 3 , 7ct2 = ei T r . (6.32) 

The kappa symmetry matrix T K for the embedding at hand is: 

e <t> 

T K = — r^o...^ 7 CTl(T2 . (6.33) 

Moreover, by using the projection conditions satisfied by the Killing spinors e, one can prove 
that 

p+q+s ± ~ p+q+s ± 



e 2 r r f 3 e = e2 ^ coso;r 1 2 + sina;r 1 f2) e , (6.34) 
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and, since the determinant of the induced metric is = p+< ^ +s , it is immediate to 

verify that the kappa symmetry projection T K e = e coincides with the projection satisfied 
by the Killing spinors of the background. Therefore, our brane probe preserves all the 
super symmetries of the background. Notice that the induced metric on the worldvolume 
along the a±, a 2 directions has the form 



2 



el[ ip + q 4 +s) do* + do*] , (6.35) 

which is conformally equivalent to the metric of a cylinder. After a simple calculation one 
can prove that the energy density of these solutions is 



H = d r 



e 2<t> 



(^pr + (q + s - p) ^ 



/coth2r 



2 sinh 2 2r 



(6.36) 



One can also have cylinders located at the south pole of the (9, ip) and (9, ip) two-spheres. 
Indeed, the above equations remain valid if 9 = n (9 = n) if one changes p — * — p (q — > —q, 
respectively). On the other hand, if p — the angle ip is constant and, as the pullback of 
the e* frame one- forms also vanishes when 9 is also constant, it follows that 9 can have any 
constant value when p = 0. Similarly, if q — one necessarily has <p — <p and 9 can be an 
arbitrary constant in this case. 

When p = 1, q = n and s = 2m, the angular part of the embedding is the same as in 
the (n, m)-winding solitons. Actually, these cylinder solutions correspond to formally taking 
r* — > — oo is the abelian solution of eq. (5.37). This forces one to take 9 = 0,7r and, thus 
one can regard the cylinder as a zoom which magnifies the region in which the probe goes 
to infinity. One can also get cylinder embeddings by consider the limit of the non-abelian 
solutions in which the probe reaches the point r = 0. For example, by taking r* = in 
eq. (6.15) one gets the p — q — 1, s — cylinder solutions while the r* — > limit of the 
embedding (6.27) corresponds to a cylinder with p — 1 and q = s = 0. Actually, when one 
takes the r* — > limit of these non-abelian embeddings one obtains two cylinder solutions, 
one with 9 = and the other with 9 = ir. This suggests that, in order to obtain a consistent 
solution, one must combine in general two cylinders located at each of the two poles of the 
(9, ip) two-sphere. Notice that this is also required if one imposes the condition of RR charge 
neutrality of the two-sphere at infinity. 



7 Quadratic fluctuations around the unit-winding em- 
bedding 

As mentioned in the introduction, we are now going to consider some of the brane probe 
configurations previously found as flavor branes, which will allow us to introduce dynamical 
quarks in the Af = 1 SYM theory. Following refs. [19, 20, 21], the spectrum of quadratic 
fluctuations of the brane probe will be interpreted as the meson spectrum of Af = 1 SQCD. 
So, let us try to elaborate on the reasons to consider these probes as the addition of flavors 
to the field theory dual. In fact, when considering the 't Hooft expansion for large number of 
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colors, the role of flavors is played by the boundaries of the Feynman graph. From a gravity 
perspective, these boundaries correspond to the addition of D-branes and open strings in the 
game. 

In our case, we have a system of N D5-branes wrapping a two-cycle inside the resolved 
conifold and Nf D5-branes that wrap another two-submanifold, thus introducing Nf flavors 
in the SU(N) gauge theory. Taking the decoupling limit with g s a'N fixed and large is 
equivalent to replacing the N D5-branes by the geometry they generate (the one studied 
in section 2), while the Nf-D5 branes that do not backreact (because we take Nf much 
smaller than N) are treated as probes. From a gauge theory perspective, this is equivalent 
to consider the dynamics of gluons and gluinos coupled to fundamentals, but neglecting the 
backreaction of the latter. Of course it would be of great interest to find the backreacted 
solution. 

The way of adding fundamental fields in this gauge theory from a string theory perspective 
was discussed in [28], where two possible ways, adding D9 branes or adding D5 branes as 
probes, were proposed. In this paper we are considering the cleaner case of D5 probes. A 
careful analysis of the open string spectrum shows the existence of a four dimensional gauge 
J\f = 1 vector multiplet and a complex scalar multiplet. This is the spectrum of SQCD. In 
the case analyzed below we will consider abelian DBI actions for the probes, so that we will 
be dealing with the Nf = 1 case. 

We have found several brane configurations in the non-abelian background which, in 
principle, could be suitable to generate the meson spectrum. One of the requirements we 
should demand to these configurations is that they must incorporate some scale parameter 
which could be used to generate the mass scale of the quarks. Within our framework such a 
mass scale is nothing but the minimal distance between the flavor brane and the origin, i.e. 
what we have denoted by r*. This requirement allows to discard the cylinder solutions we 
have found since they reach the origin and have no such a mass scale. We are thus left with 
the unit-winding solutions of section 6.1 and the zero-winding solutions of section 6.2 as the 
only analytical solutions we have found for the non-abelian background. 

In this section we shall analyze the fluctuations around the unit-winding solutions of sec- 
tion 6.1. We have several reasons for this election. First of all, the unperturbed unit- winding 
embedding is simpler. Secondly, we will show in appendix B that the UV behaviour of the 
fluctuations is better in the unit-winding configuration than in the zero-winding embed- 
ding. Thirdly, the unit-winding embeddings of constant ip incorporate the correct pattern 
U(l) — > 7L 2 of R-symmetry breaking, whereas for the zero-winding embeddings of eq. (6.27) 
the U(l) symmetry is broken to a S 4 subgroup. 

Recall from section 6.1 that we have two possible solutions with ip = 0,ir (mod27r), 
which are the ones displayed in eqs. (6.15) and (6.16). As discussed in section 6.1, the 
solution of eq. (6.16) contains a one-parameter subfamily of embeddings which reach the 
origin and, thus, they should correspond to massless dynamical quarks. On the contrary, 
the embeddings of eq. (6.15) pass through the origin only in one case, i.e. when r* = 
and, somehow, the limit in which the quarks are massless is uniquely defined. Recall that 
for r* = the solution (6.15) is identical to the unit-winding cylinder. For these reasons we 
consider the configuration displayed in eq. (6.15) more adequate for our purposes and we 
will use it as the unperturbed flavor brane. 
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We will consider first in subsection 7.1 the fluctuations of the scalar transverse to the 
transverse probe, while in subsection 7.2 we will study the fluctuations of the worldvolume 
gauge field. The gauge theory interpretation of the results will be discussed in subsection 
7.3 



7.1 Scalar mesons 

Let us consider a non-abelian unit-winding embedding with 9 — 9, <p — tp + constant and 
■ip = 7r (mod27r). For convenience we take first r as worldvolume coordinate and consider 9 
as a function of r, ip and of the unwrapped coordinates x, i.e. 9 = 9(r, x, ip). The lagrangian 
density for such embedding can be easily obtained by computing the induced metric. One 
gets: 

£ = -e 2 * sin# x 



x 



l + rtanhr((<9 r #) 2 + {d x 9f 



+ 



1 



cos 2 9 + r coth r sin 9 



+ rd r 9 - cot 9 



(r cothr + cot 2 9) + 



(7.1) 



where we have neglected the term d r (r cos 9e 2( ^) which, being a total radial derivative, does 
not contribute to the equations of motion. 

We are going to expand this lagrangian around the corresponding non-abelian unit- 
winding configuration obtained in section 6.1. Actually, by taking in eq. (6.15) rj — +1 one 
obtains a configuration with tp = n (mod27r), which corresponds to a function 9 = 9 (r), 
given by: 

sin^ (r) = ^, (7.2) 
smhr 

where r* is the minimum value of r and r* < r < oo. It is clear from this equation that with 
the coordinate r we can only describe one-half of the brane probe: the one that is wrapped, 
say, on the north hemisphere of the two-sphere, in which 9 e (0,7r/2). Outside this interval 
^o(^) is a double-valued function of r. Let us put 



9(r,x,ip) = 9 (r) + x(r,x,<p) , 



(7.3) 



and expand C up to quadratic order in x- Using the first-order equation satisfied by 6*o(r), 
namely: 

d r 9 = — cothr tan 9 , 



(7.4) 



we get 
£ = 



2 1 + r coth r tan 2 9 
1 



i n /a \2 2r r cothr 2 

r t&nh r cos 9 o (OrX) + - x &rX + 



cos 9 



cos 3 9 



X 



— — e 2 ^ r tanh r 



cos9 (d x x) + 



cos # ( 1 + r coth r tan 2 9 ) 



(d v xf 



(7.5) 
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In the equations of motion derived from this lagrangian we will perform the ansatz: 

X (r,x,<p) = e ikx e il * £(r) , (7.6) 

where, as tp is a periodically identified coordinate, / must be an integer and k is a four-vector 
whose square determines the four-dimensional mass M of the fluctuation mode: 

M 2 = -k 2 . (7.7) 

By substituting the functions (7.6) in the equation of motion that follows from the lagrangian 
(7.5), one gets a second order differential equation which is rather complicated and that can 
only be solved numerically. However, it is not difficult to obtain analytically the asymptotic 
behaviour of £(r). This has been done in appendix B and we will now use these results to 
explore the nature of the fluctuations. For large r, i.e. in the UV, one gets (see eq. (B.12)) 
that £(r) vanishes exponentially in the form: 

f (r) ~ e — cos [VM 2 -P r + 5} , (r -> oo) , (7.8) 

f 4 

where 5 is a phase and we are assuming that M 2 > I 2 . For M 2 < I 2 the fluctuations do not 
oscillate in the UV and we will not be able to impose the appropriate boundary conditions 
(see below). Notice that our unperturbed solution 9 (r) also decreases in the UV as 

6 (r) ~ e~ r (r -> oo) . (7.9) 

Thus C( r )/^o( r ) ~~ * as r — > oo and the first-order expansion we are performing continues 2 to 
be valid in the UV. On the other hand, for r close to r* there are two independent solutions, 
one of them is finite at r = r* while the other diverges as 

«r) ~ ^JU . (7.10) 

Let us now see how one can use the information on the asymptotic behaviour of the 
fluctuation modes to extract their value for the full range of the radial coordinate. First of 
all, it is clear that, in principle, by consistency with the type of expansion we are adopting, 
one should require that £ << 6q. Thus, one should discard the solutions which diverge in the 
infrared (see, however, the discussion below). Moreover, the behaviour of the fluctuations 
£ for large r should be determined by some normalizability conditions. The corresponding 
norm would be an expression of the form: 



/ dr^ e , (7.11) 
j o 



where v /7 is some measure, which can be determined by looking at the lagrangian (7.5). If 
we regard x as a scalar field with the standard normalization in a curved space, then is 
just the coefficient of the kinetic term | (d r x) 2 in £, namely 



2( p rtanhrcos^o 



v 1 + r coth r tan 2 6»o ' 

2 For the zero-winding solution, on the contrary, the ratio £{r)/6o(r) diverges in the UV (see appendix B). 
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For large r, ^7 behaves as 



( r 



00 ) 




(7.13) 



Notice that the factors on the right-hand side of eq. (7.13) cancel against the exponentials 
and power factors of £ 2 in the UV (see eq. (7.8)). As a consequence, all solutions have 
infinite norm. 

The reason for the bad behaviour we have just discovered is the exponential blow up of 
the dilaton in the UV which invalidates the supergravity approximation. Actually, if one 
wishes to push the theory to the UV one has to perform an S-duality, which basically changes 
e 2 ^ by e~ 2< ^. The S-dual theory corresponds to wrapped Neveu-Schwarz fivebranes and is 
the supergravity dual of a little string theory. Notice that, by changing e 2 ^ e~ 2< ^ in the 
measure (7.12), all solutions become normalizable, which is as bad as having no normalizable 
solutions at all. Moreover, it is unclear how to perform an S-duality in our D5-brane probe 
and convert it in a Neveu-Schwarz fivebrane for large values of the radial coordinate. 

A problem similar to the one we are facing here appeared in ref. [30] in the calculation 
of the glueball spectrum for this background. It was argued in this reference that, in order 
to have a discrete spectrum, one has to introduce a cut-off to discriminate between the two 
regimes of the theory. Notice that, since they extend infinitely in the radial direction, we 
cannot avoid that our D5-brane probe explores the deep UV region. However, what we can 
do is to consider fluctuations that are significantly nonzero only on scales in which one can 
safely trust the supergravity approximation. In ref. [30] it was proposed to implement this 
condition by requiring the fluctuation to vanish at some conveniently chosen UV cut-off A. 
Translated to our situation, this proposal amounts to requiring: 



This condition, together with the regularity of £(r) at r = r*, produces a discrete spectrum 
which we shall explore below. Notice that, for consistency with the general picture described 
above, in addition to having a node at r = A as in eq. (7.14), the function £(r) should be 
small for r close to the UV cut-off. This condition can be fulfilled by adjusting appropriately 
the mass scale of our solution, i. e. the minimal distance r* between the probe and the origin, 
in such a way that r* is not too close to A. 

Notice that, by imposing the boundary condition (7.14) on the fluctuations, we are effec- 
tively introducing an infinite wall located at the UV cut-off. The introduction of this wall 
allows to have a discrete spectrum and should be regarded as a physical condition which 
implements the correct range of validity of the background geometry as a supergravity dual 
of M = 1 Yang-Mills. Even if this regularisation could appear too rude and unnatural, the 
results obtained by using it for the first glueball masses are not too bad [30] . 

The cut-off scale A should not be a new scale but instead it should be obtainable from 
the background geometry itself. The proposal of ref. [30] is to take A as the scale of 
gaugino condensation, which is believed to correspond to the point at which the function 
a(r) approaches its asymptotic value a = 0. A more pragmatic point of view, to which we 
will adhere here, is just taking the value of A which gives reasonable values for the glueball 
masses. In ref. [30] the value A = 2 was needed to fit the glueball masses obtained from 
lattice calculations, whereas with A = 3.5 one gets a glueball spectrum which resembles 




(7.14) 
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Figure 4: Graphic representation of the first three fluctuation modes for r* = 0.3, A = 3 and 
/ = 0. The three curves have been normalized to have £(r*) = 1. 

that predicted by other supergravity models. Notice that from r = to r = 3 the effective 
string coupling constant increases in an order of magnitude. From our point of view it is 
also natural to look at the effect of the background on our brane probe. In this sense it is 
interesting to point out that for r* = 2 — 3 onwards the abelian and non-abelian embeddings 
are indistinguishable (see sect. 6.1). 

We have performed the numerical integration of the equation of motion of £(r) subject 
to the boundary condition (7.14) by means of the shooting technique. For a generic value 
of the mass M the solution diverges at r = r*. Only for some discrete set of masses the 
fluctuations are regular in the IR. In figure 4 we show the first three modes obtained by this 
procedure for = 0.3 and A = 3. From this figure, we notice that the number of zeroes 
of £(r) grows with the mass. In general one observes that the n th mode has n — 1 nodes in 
the region r* < r < A, in agreement with the general expectation for this type of boundary 
value problems. Moreover, for / = 0, the mass M grows linearly with the number of nodes 
(see below). 

At this point it is interesting to pause a while and discuss the suitability of our election 
of r as worldvolume coordinate. Although this coordinate is certainly very useful to extract 
the asymptotic behavior of the fluctuations (specially in the UV), we should keep in mind 
that we are only describing one half of the brane, i.e. the one corresponding to one of the two 
hemispheres of the two-sphere. On the other hand, the election of the angle as excited scalar 
has some subtleties which we now discuss. Actually, to describe the displacement of the 
brane probe with respect to its unperturbed configuration it is physically more sensible to 
use the coordinate y, defined in eq. (5.20). Accordingly, let us define the function y(r,x,<p) 
as 

y(r,x,ip) = r sin 9(r, x, ip) . (7-15) 

Let us put in this equation 9(r,x,ip) = 0o(r) + xij'^iV)- At the linear order in x we are 
working, y can be written as 

y(r,x,<p) = yoij-) + r cos 6 (r)x(r,x, if) , (7.16) 

where yo(r) = r sin 6> (r) corresponds to the unperturbed brane. Notice, first of all, that for 
those modes in which x( r *, x , f) is finite, the fluctuation term in y(r*,x,(p) vanishes since 
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cos 6*o (r) — > when r — > r*. Then 

y(r*,x,ip) = y (n) = r* if <£>) is fi nite • ( 7 -17) 

Thus, by considering those modes x that are regular at r = r* we are effectively restricting 
ourselves to the modes which have a node in the y coordinate at r = r*. If, on the contrary 
X diverges for r « r*, we know from eq. (7.10) that it behaves as x ~ W r — r *- But we 
also know that cos#o(r) - ► when r — > r* and, in fact, the second term on the right-hand 
side of eq. (7.16) remains undetermined. The precise form in which cos#o( r *) vanishes can 
be read from eq. (B.13), namely cos 9 (r) — r* for r k r». Therefore, even if x 

diverges at r = r*, we could have, in the linearized approximation we are adopting, a finite 
value for y(r*,x,tp). Thus, modes with x( r *> x i V 7 ) ~~ * 00 should not be discarded. Actually, 
y(r*,x,(p) — yo( r *)) although finite, is undetermined in eq. (7.16) for these modes and, in 
order to obtain its allowed values we should impose a boundary condition at the other half 
of the brane. As previously mentioned, this cannot be done when r is taken as worldvolume 
coordinate. Therefore, it is convenient to come back to the formalism of sects. 3-5, in which 
had been chosen as one of the worldvolume coordinates. In this approach, the unperturbed 
brane configuration is described by a function r (9), given by: 

sinhro(tf) = ^ , (7.18) 

and the brane embedding is characterized by a function r = r(9,x,f), which we expand 
around r (9) as follows: 

r(9,x,<p) = r (9) + p(6,x,<p) . (7.19) 

Plugging this expansion into the DBI lagrangian of eq. (5.24) and keeping up to second 
order terms in p, one gets the following lagrangian density: 
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7.20) 



Similarly to what we have done with the lagrangian (7.5), we will look for solutions of the 
equations of motion of C which have the form: 

p(6,x,<p) = e^e^CW , (7.21) 

where / is an integer and k 2 = —M 2 . As before, in order to get a discrete spectrum one must 
impose some boundary conditions. In the present approach we should cutoff the regions 
close to the two poles of the two-sphere. Accordingly, let us define the following angle 
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Figure 5: Plot of 5y = £sin# versus x = rcos9 for the first four modes for r* = 0.3, A = 3 
and / = 0. The dashed and dotted curves pass through the origin and correspond to the first 
two modes of figure 3. 

where, as indicated, we are taking 9\ in the range < 9 A < vr/2. Notice that 9 A — > if 
A — > oo as it should. Clearly 9\ and ix — 9\ are the two angles that correspond to the radial 
scale A. Therefore, we impose the following boundary conditions to our fluctuation 

C(0a) = C(vr - 0a) = . (7.23) 

The equations of motion derived from (7.20), subjected to the boundary conditions (7.23) 
can be integrated numerically by means of the shooting technique. One first enforces the 
condition at 9 = 9^ and then varies the mass M until ((n — 9^) vanishes. This only happens 
for a discrete set of values of the mass M. For a given value of I, let us order the solutions 
in increasing value of the mass. In general one notices that the n th mode has n — 1 nodes 
in the interval 9\ < 9 < n — 9^ and for n even (odd) the function ((9) is odd (even) under 
9 — > 7r — 9. In figure 5 we have plotted the first four modes corresponding to r* = 0.3, A = 3 
and I = 0. The modes odd under 9 — > n — 9 vanish at 9 = it/2 and their masses and shapes 
match those found with the lagrangian (7.5) and the boundary condition (7.14). On the 
contrary, the modes with an even number of nodes in 9\ < 9 < ix — 9\ are the ones we were 
missing in the formulation in which r is taken as worldvolume coordinate. 

Let M n j(r*,A) be the mass corresponding to the n th mode for a given value / and the 
mass scales r* and A. Our numerical results are compatible with an expression of M n ^(r*, A) 
of the form: 

M n ,,(r*,A) = y/m?(r„A) n 2 + P (7.24) 

To illustrate this fact we have plotted in figure 6 the values of M n ^(r*,A) for r* = 0.3 and 
A = 3, together with the curves corresponding to the right-hand side of (7.24). 

We have also studied the dependence of the coefficient m(r*, A) on the two scales (r*, A). 
Recall that r* is the minimal separation of the brane probe from the origin and, thus, can 
be naturally identified with the mass of the quarks. We obtained that ra(r*,A) can be 
represented by an expression of the type: 

m(r„A) = ^ + KAK 2 . (7.25) 

The coefficients on the right-hand side of eq. (7.25) have been obtained by a fit of m(r*, A) 
to a quadratic expression in r*. An example of this fit is presented in figure 7. The first 
term in eq. (7.25) is a universal term (independent of r*) which can be regarded as a 
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Figure 6: Mass spectrum for r* = 0.3 and A = 3 for I — 0, • • • , 6. The solid lines correspond 
to the right-hand side of eq. (7.24). The triangles (squares) are the masses of the modes 
((9) which are even (odd) under 9 — > n — 9. 
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finite size effect induced by our regularisation procedure. We also have determined the 
dependence of the coefficient b on A and it turns out that one can fit 6(A) to the expression 
6(A) = 0.23 A" 2 + 0.53 A" 3 

We have obtained a very regular mass spectrum of particles, classified by two quantum 
numbers n, I (see eqs. (7.24) and (7.25)). We can offer an interpretation of these formulas. 
Indeed, recall that the two-submanifold on which we are wrapping the brane is topologically 
like a cylinder, with the compact direction parametrized by the angle tp. The quantum 
number / is precisely the eigenvalue of the operator —id^, which generates the shifts of tp 
and, indeed, the dependence on I of the mass displayed in eq. (7.24) is the typical one 
for a Kaluza-Klein reduction along a compact coordinate. Therefore, we should interpret 
the mesons with I ^ as composed by Kaluza-Klein modes. Moreover, the term in (7.25) 
proportional to r 2 can be understood as the contribution coming from the mass of the 
constituent quarks, while the term ^ can be interpreted as a contribution coming from the 
'finite size' of the meson. Indeed, it looks like a Casimir energy and it is originated in the 
presence of the cut-off A. 

It is perhaps convenient to emphasize again that the cut-off procedure implemented here 
is a very natural procedure for this type of computations in this supergravity dual. In fact, 
given that the mesons are an IR effect in SQCD, we expect the contributions of high energy 
effects to be irrelevant or negligible in the physical properties of the meson itself. This 
is indeed what we are doing when cutting off the integration range. We are just taking 
into account the 'non-abelian' part of the probes, while neglecting the 'abelian' part or, 
equivalently, discarding high energy contributions. 



7.2 Vector mesons 



Let us now excite a gauge field that is living in the worldvolume of the brane. The linearized 
equations of motion are: 



On 



V^tf nm = o , 



(7.26) 



where g st is the determinant of the induced metric in the string frame and T nm is the field 
strength of the worldvolume gauge field A m , i.e. T nm = d n A m — d m A n . Let us assume 
that the only non-vanishing components of the gauge field A are those along the unwrapped 
directions of the worldvolume of the brane. In what follows we are going to use 9 as 
worldvolume coordinate. Let us put 



(7.27) 



where is a constant polarization four-vector and I must be an integer. It follows from the 
equations of motion (7.26) that must be transverse, i.e.: 



fc% = o , 



and that q(9) must satisfy the following second-order differential equation: 
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where, as in eq. (7.7), k 2 = —M 2 . We have solved numerically this equation by means of 
the shooting technique with the boundary conditions 

?(0a) = <^ - A ) = . (7.30) 

Surprisingly, the set of possible values of M is given by the same expression as in the scalar 
meson case, i.e. eq. (7.24), with a coefficient function m(r*, A) which is, within the accuracy 
of our approximate calculation, equal numerically to that of the scalar mesons. This is quite 
remarkable since the differential equations we are solving in both cases are quite different (the 
equation satisfied by ((9) is much more complicated than eq. (7.29)). Thus, to summarize, 
we predict a degeneracy between the scalar and vector mesons in the corresponding Af = 1 
SYM theory. 

7.3 Gauge theory interpretation 

Let us comment on some gauge theory aspects that can be read from the brane probes in 
sections 3 to 7. For this purpose, we shall concentrate on the main solutions that have 
been used, namely, the abelian solution with unit winding (n = 1) in eq. (5.22) and its 
non-abelian extension of eq. (6.15). 

First of all, let us analyze the R-symmetry of the gauge theory from the probe viewpoint. 
It is clear that the abelian solution (5.22) is invariant under shifts of ip by a constant since 
the value of this angle is an arbitrary constant in this solution. This symmetry has been 
identified as the geometric dual of the R-symmetry [6, 34, 35]. Actually, this is not a U(l) 
invariance of the background because [6, 34, 35] of the presence of the RR form that selects, 
by consistency, only some particular values of the angle ip, i.e. ip = ^ with 1 < n < 2N. So, 
the abelian probes can see a TL 2N symmetry. In contrast, when we consider the non-abelian 
probe, the solution (6.15) selects two particular values of the angle ip, thus breaking 7L 2 n 
down to TL 2 . Notice that the U(l) — > TL 2 n breaking is an UV effect (it takes place already in 
the abelian background) while the TL 2 n —* TL 2 breaking is an IR effect which appears only 
when one considers the full non-abelian regular background. This same breaking pattern was 
observed in the case of SQCD with massive flavors. Indeed, as showed in [29], the theory 
with massive flavors has a non-anomalous discrete TL 2N R-symmetry, given by (component 
fields are used, A is the gaugino and $ and $ are the squarks): 

A -> e - im/N \ , $ -> e - i7Tn/N <f> , $ -> e -' mn/N $ , (7.31) 

with n = 1,...2N 

As shown in [29], this TL 2N symmetry is broken down to TL 2 by the formation of a squark 
condensate. Indeed, one can see that < > transforms as < >^ e -*™/N < > 
leaving us with a TL 2 . Besides, the squark condensate is consistent with supersymmetry, 
because the F-term equation of motion < > — m = is satisfied. Notice that this 
preservation of symmetry when m ^ is in agreement with the kappa symmetry of our 
brane probes. 

Apart from all this, there is a vectorial U(l) symmetry that remains unbroken in our 
brane probe analysis and that we have associated with the invariance under translations in 
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(p. On the field theory side this symmetry can be identified with a phase change of the full 
chiral/antichiral multiplet $ — > e ta §, $ — > e~*°5>, which is nothing but the U(1)b baryonic 
number symmetry of the theory. The two possible assignations ±1 of the baryonic charge 
are in correspondence with the two possible identifications of the S 1 described by (p with the 
S 1 parametrized by if. The spectrum we have found is independent of this identification. 

We would like also to comment briefly on the possibility of taking the parameter r* = 0. 
Given that this parameter can be associated with the mass of the quark (since it is the 
characteristic distance between the probe brane and the background) one would like to study 
the case in which this parameter is taken to be zero. Nevertheless the approach implemented 
here seems to break down for this particular value of the parameter. Indeed, taking r* = 
will imply that 9 = 0, so, a fluctuation of it can lead to negative values of 9 taking us 
out of the range of this coordinate. The fact that our approach apparently does not work 
for the case of massless quarks seems to be in agreement with the fact that SQCD with 
massless flavors, has some special properties like spontaneous breaking of supersymmetry, 
the existence of a runaway potential (Affleck-Dine-Seiberg potential) and the non-existence 
of a well-defined vacuum state. Notice that, when r* ^ 0, our approach, by construction, 
deals with massive quarks that preserve supersymmetry, since our probes were constructed 
by that requirement. 

8 Summary and conclusions 

Let us summarize and repeat here some points that we believe are worth stressing again. 

We studied a supergravity dual to M = 1 SYM, based on a geometry representing a 
stack of D5 branes wrapping a two-cycle. In this gravity solution, we have found special 
surfaces where one can put probe D5-branes without spoiling supersymmetry. This was 
done, basically, by imposing that the Killing spinors of the background satisfy the projection 
imposed by the kappa symmetry of the D5-brane probe. 

A wide variety of kappa symmetric solutions was found and a rich mathematical struc- 
ture was pointed out. In fact, solutions where the probes are at a fixed distance from the 
" background branes" are shown to break supersymmetry. This phenomenon is in agreement 
with the non-existence of moduli space of Af = 1 SYM. On the other hand, we studied solu- 
tions corresponding to the abelian background (that is to say the large r regime of the full 
background). These abelian solutions are shown to have a very interesting analytic structure 
(harmonic functions and Cauchy-Riemann equations show up) allowing interesting explicit 
solutions. In extending these studies to the full background, we found several classes of non- 
abelian solitons. We think that it might be possible that these solutions show themselves 
useful when studying other aspects of the model not explored in this paper. 

Then, in section 7, we have used the kappa symmetric solutions (in particular those that 
we called unit-winding) to introduce fundamental matter in the dual to SYM theory. 

Given that the brane probes do not backreact on the background, these flavors are intro- 
duced in the so called quenched approximation. Nevertheless, many qualitative features and 
quantitative predictions for the strong coupling regime of M = 1 SQCD can be addressed. 
Among them, the qualitative difference between the massless and massive flavors are clear 
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in this picture. Indeed, by construction, our approach deals with the massive-flavor case, 
so, the problems or peculiarities of the massless case can be seen in this approach under a 
different geometrical perspective. 

Other characteristic feature of SQCD with few flavors is the breaking pattern of R- 
symmetry. The 7L 2N —> ^2 breaking is geometrically very clear from the brane probe 
perspective. Also, the preserved U(1)b is geometrically realized, as explained in the previous 
section, by arbitrary changes in the coordinate (p. 

A mass spectrum for the low energy excitations (mesons) of massive SQCD was found 
and we believe this is a very interesting and quantitative prediction of this paper. In fact, a 
nice formula for the masses is derived that exhibits a BPS-like behavior with the level (n) 
and with the Kaluza-Klein quantum number (/) of the meson. Basically, our formula gives 
the meson masses in terms of the mass of the fundamentals and the Casimir energy due to 
finite size of the meson, and shows explicitly the contamination of the meson spectrum due 
to the Kaluza-Klein modes. It would be very interesting if lattice calculations could validate 
or invalidate the formula found here. 

The appendices provide a hopefully very useful technical information for other workers 
in the area which contain, in particular, the explicit expressions of the Killing spinors of the 
background considered here and for another related background dual to M = 1 SYM theory 
with a Chern-Simons term in 2 + 1 dimensions. 

We believe that this paper opens some interesting routes of future research. To mention 
some of them, it should be interesting to repeat the procedure here in other "wrapped 
branes" setups, dual to lower supersymmetric field theories, to see if the features for the 
mass formula and symmetries are preserved. In particular, it would be very interesting to 
apply the methods employed here to the Klebanov-Strassler background [23]. The analysis of 
ref. [36] of the glueball spectrum suggests that, contrary to what happens to the background 
studied here, the spectrum of mesons could be obtained without introducing any additional 
cut-off. Actually, this meson spectrum was obtained in ref. [22] from the fluctuations of a 
D7-brane probe wrapped on a submanifold which is analogous to our cylinders. Presumably, 
our method would allow us to obtain a family of supersymmetric embeddings which depend 
on a parameter that can be identified with the distance between the brane probe and the 
origin. 

On the other hand, the brane probe lagrangian should be regarded as a sort of "low energy 
effective lagrangian" for the strong coupling limit of SQCD with few flavors; it would be very 
nice to develop the physics of this effective theory because, possibly, the scattering of mesons 
and other interesting features, like low energy theorems, could be studied or understood in 
this geometrical context. 

Extending the treatment of this paper to non supersymmetric theories might seem a 
contradiction (the paper is based on kappa symmetry) but it is likely that with softly broken 
backgrounds like those appeared recently in the literature, see ref. [12], similar methods to 
the ones exhibited here can be used. 
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A Supersymmetry of the J\f = 1 supergravity duals 

The supersymmetry transformations for the type IIB dilatino A and gravitino \& M when the 
RR three-form is nonzero are: 

The Killing spinors of the background are those e for which the right-hand side of eq. (A.l) 
vanishes. In order to satisfy the equations 5X = = we will have to impose certain 
projection conditions to e. Once these projections are imposed, the equations 5\ = 5^/^ = 
are equivalent to a system of first-order differential equations for the metric and RR three- 
form. Moreover, following the methodology of ref. [16] (see also ref. [8]), we will be able to 
determine the explicit form of e from the projections that satisfies. In this appendix we will 
carry out this program for two Af = 1 supergravity duals which correspond to D5-branes 
wrapping a two- and a three-cycle. 

A.l Fivebranes wrapped on a two-cycle 

For a metric ansatz such as that of eq. (2.1), let us consider the frame 

e xi = e idx i , (i = 0,1,2,3) , 

e 1 = e i +h d0, e 2 = ei +h sin 9dip , 

e r = e%dr, e l = ^ ( w i - A 1 ) , (i = 1,2,3). (A.2) 
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In order to solve the equations 5 X = = 0, we impose to the spinor e the condition for a 
SUSY two-cycle: 

ri 2 e = f 12 e, (A.3) 

and the following projection 

e = ie* . (A.4) 
Then, the condition 5X = yields the equation 

e~ 2h . 1 

<t>'e - (1 + —(a 2 -l))T r T 123 e + -a' e~ h Y 2 Y 2 e = . (A.5) 

Let us now consider the supersymmetric variation of the different components of the grav- 
itino. It is easy to verify that, with the projections (A.3) and (A.4), the equation corre- 
sponding to the unwrapped directions x l and that corresponding to the directions i, coincide 
with the one obtained from the variation of the dilatino (eq. (A.5)). On the other hand the 
variation of the components along the two-sphere gives rise to the following equation 

tie - i a' e- h T 2 t 2 e - a e~ h r 2 f 2 r r f 123 e + l - (a 2 - 1) e~ 2h T r t 123 e = , (A.6) 

whereas, after using eq. (A.5), the condition 5tp r = is equivalent to 

d r e - \a'e- h V 2 t 2 e - \<p' e = . (A.7) 
2 8 

Moreover, from the variation of the dilatino (eq. (A.5)) we learn that e must satisfy 

r r f 123 e = ((3 + pr 2 r 2 )e, (A.8) 

where (3 and (5 can be read from eq. (A.5), namely 

(f)' ~ _ 1 e"V 

1 + ^ (a 2 - 1) ' 2l + ^(a 2 -l) ' (,) 

On the other hand, from the condition (r r fi 23 ) 2 e = e and since {r r fi 23 , T 2 f 2 } = 0, it follows 
that: 

[3 2 + f3 2 = l, (A.10) 
and, therefore, we can represent (5 and (3 as: 

P = cosa, $ — sma. (A. 11) 

Substituting the radial projection (eq. (A.8)) into eq. (A.6) and considering the terms with 
and without T 2 r 2 , we get the following two equations: 

ti = --e- 2h {a 2 - l)(3-ae~ h (3 , a' = -2a/3 + e~ h {a 2 - 1) (3 . (A.12) 
2 



40 



By using the definition of (3 and (3 (eq. (A. 9)) into the second equation in (A. 12), we get the 
following relation between 0' and a': 



a 



0' = - — 
r 2a 



1) 



(A.13) 



Furthermore, from the condition (3 2 + (3 2 = 1 one obtains a new relation between 0' and a', 
namely: 

/2 + ^e~ 2h a' 2 = [1 + I e -^(a 2 -l)] 2 . (A.14) 

By combining eqs. (A.13) and (A.14) one can get the expression of 0' and a' in terms of a 
and h. Moreover, by using these results in eq. (A. 9), one can get (3 and (3 as functions of 
a and h and, by plugging the corresponding expressions on the first eq. in (A. 12), one can 
obtain the differential equation for h. In order to write these expressions in a compact form, 
let us define: 



Q 



l e 4h + La(„2 + 1 ) + _L( a 2_!)2 

2 16 



(A.15) 



Then, one has the following system of first-order differential equations for 0, h and a: 



ti = 



a = 



1 

Q 
l 

2Q 
2a 



-2h 



e 2h _ 



16 



(a 2 - If 



-2h 



a 2 + l + 



[a 2 - If 



e 2h + -Aa 2 - 1) 
4 



and the values of (3 — cos a and (3 = sin a, which are given by: 



sin a = 



ae 



Q 



J2h 



cos a = 



Q 



(A.16) 



(A.17) 



It is interesting to notice that, when solving the quadratic eq. (A.14) to obtain (A.16) and 
(A.17), we have a sign ambiguity. We have fixed this sign by requiring that h! is always 
positive. 

It remains to verify the fulfillment of equation (A. 7). Notice first of all that the radial 
projection (A. 8) can be written as 



i r 1 123 e — e e 



which, after taking into account that {T r f 123, T 2 f 2} = 0, can be solved as: 

Tr Tl23 e o — e o • 



(A.18) 



(A.19) 



Inserting this parametrization of e into eq. (A. 7), we get two types of terms, with and 
without r 2 r2 ? which yield the following two equations: 



d r e - 



, 



a 



-e~ h a' 



(A.20) 
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The first of these equations can be solved immediately, namely 

eo = e%7] , (A.21) 



with i] being a constant spinor. Moreover, by differentiating eq. (A. 17) and using eq. (A. 16), 
one can verify that the second equation in (A. 20) is satisfied. Thus, the explicit form of the 
Killing spinor is 

e = e t Flfl e* 7] , (A.22) 
where r\ is a constant spinor satisfying: 

T x o... x iT l2 ri = r} , T 12 r] = f 12 r] , r} = irf . (A.23) 

Eq. (A. 16) is a system of first-order differential equations whose solution determines the 
metric, dilaton and RR three-form of the background. It can be verified that the values 
given in eq. (2.4) solve the system (A. 16). Moreover, by plugging in eq. (A. 15) the values 
of h and a given in eq. (2.4), one verifies that 

Q = r . (A.24) 

Then, one gets the following simple expression for cos a: 

2r 

cos a = coth2r — 5 — . (A. 25) 

sinh 2 2r v ; 

On the other hand, being a spinor of definite chirality of type IIB supergravity, e satis- 
fies T x o... x iTi2T r Ti23e = e. Thus, if we multiply the radial projection condition (A. 8) by 
T x o... x aTi2, we obtain: 

r x o... x 3 ( coso;ri2 + sinaTir^e = e . (A. 26) 

Eqs. (A. 3), (A. 4) and (A. 26) are the ones needed in the analysis of the kappa symmetry of 
a D5-brane probe (eq. (2.12)). 

A. 2 Fivebranes wrapped on a three-cycle 

As a second example of the calculation of the Killing spinors of a Af = 1 supergravity dual, 
let us consider the supergravity solution for a D5 brane wrapped on a three-sphere. In this 
case the ansatz for the metric in Einstein frame is 

ds 2 = e 2 [ dx\ 2 + ^ R(r) 2 ( a* f + dr 2 + ^ ( w l - A 1 f } , (A.27) 

where is the dilaton, R(r) is a function to be determined and a 1 and w l (i = 1, 2, 3) are 
two sets of su(2) left invariant one-forms. The gauge field components A 1 are parametrized 
in terms of a new function uj{r) as: 

A* = 1 + " (r) a* . (A.28) 
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Let F % be the field strength of A 1 (defined as in eq. (2.6)). Its components are: 

/ 2 i 

pi = ^ dr A (j* + ^— - e ijfe ^' A a fe . (A.29) 
2, o 

The RR three-form of the background is now 

F (3) = -1 (w 1 - A 1 ) A (u> 2 - A 2 ) A (u> 3 - A 3 ) + 1 ^ F* A (iu* - A*) + /i , (A.30) 

i 

where h is determined by requiring that F( 3 ) satisfies the Bianchi identity dF^ = 0. One 
easily verifies that h must satisfy the following equation 

dh = | Yl pl A pi ■ ( A - 31 ) 

i 

By explicit calculation, one gets: 

F l A F l = - {u 2 - 1 )J e jkm dr A a j A a k A a m . (A.32) 

i 8 

Thus, the equation for h can be solved as: 

h = ^KW^^^A/ , (A.33) 

where 

V(r) = 2lu 3 - 6cu + 8A; , (A.34) 

with k being a constant. 

Let us study the supersymmetry preserved by this ansatz in the frame: 

e x * = e% dx l , (i — 0, 1, 2) , e r = e* dr , 

e* = itV\ e J = ^ (w< - A') , (i = 1,2,3), 

(A.35) 

We now impose the following projections on the spinors 

Tifie = r 2 f 2 e = r 3 f 3 e , 

e = ie* . (A.36) 
Then, the condition 5\ = becomes 

0'e " (l + 3 ^4i^) rrfW + ^ Flfie " 8^ rifirrfl236 = ° • (A ' 37) 
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Moreover, from Stpi = 0, one gets 

R'e - y r^e + T r f 123 e + ( ^ - cu ) 1^ T r f 123 e = . (A.38) 

The vanishing of the supersymmetric variation of the radial component of the gravitino gives 
rise to: 

a r e-—r 1 t 1 e--<l/e = 0, (A.39) 
where we have used eq. (A. 37). Let us solve these equations by taking the projection 

r r f 123 e = (/3 + prj^e . (A.40) 

As in case studied in section (A.l), from (r r f 12 3) 2 e = e and {T r f 12 3, TiT i} = it follows 
that /5 2 + /5 2 = 1 and thus we can take (3 = cos a and (3 = sin a. 

Let us substitute our ansatz for r r ri 23 e (eq. (A.40)) on the equations coming from the 
dilatino and gravitino (eqs. (A. 37) and (A.38)). From the terms containing the unit matrix, 
we obtain equations for 0' and Rl: 

Moreover, by considering the terms with Tif i, we obtain two expressions for u/ 

3u/ = (4i? + 3^)/3, 

a/ = (^-2w)/? + 2^0. (A.42) 

By combining these last two equations we get 

/ V \ ^ / 1 - u> 2 2i? \ ~ . . 
o -«;)/?= ( + — ) /3 . A.43 

By plugging this last relation in the condition [3 2 + (3 2 = 1 one can easily obtain the 
expression of (3 and (3. Indeed, let us define M as follows: 

f V \ 2 / 1 - W2 2R \ 2 /A ,,S 

M -(w" ra ) + (^7T + t) ■ (A ' 44) 

In terms of this new quantity M, the coefficients (3 and /3 are given by: 

1 /2R 1-lo 2 \ ~ 1 / V \ , k , 

^ = cosa = ~m { t + -w J ' ^ = sina = 7^ (24^ " "J • ( 5) 
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By using these values of (5 and $ in the equations which determine R', uj' and <fi', we obtain 
a system of first-order BPS equations which are identical to those written in refs. [31, 32, 33] 
(see also ref. [37]). They are: 



R' = 



1 



3VM 
4R 



+ ^( 3 ( 1 -^ 2 ) 2 -Vu)+u 2 + 2 



UJ 



3VM 



6AR 4 2R 2 
V 



32R 4 



1 - u z + 



2 k - uj 3 
2R 2 



— UJ 



3 M pVx 3VM 



The radial projection condition (A. 40) can be written as 

r r f 123 e = e° rifl e . 
Since {T r f 12 3, T^Y i} = 0, this equation can be solved as: 

€ = e 2 TlTl €q , T r r 12 3 eo = e . 



(A.46) 



(A.47) 



(A.48) 



Plugging now this parametrization of e into the equation obtained from the variation of the 
radial component of the gravitino (eq. (A. 39) ), we arrive at the following two equations: 



_ ~8~ C ° ' 



a = 



JR . 



The equation for eo can be solved immediately: 

e = e» T) , 



(A.49) 



(A.50) 



where 77 is a constant spinor. Moreover, one can verify that the equation for a is a consequence 
of the first-order BPS equations (A.46). Therefore, the Killing spinors for this geometry are 
of the form: 

e = e ~t r i f i e f v ; (A.51) 
where r\ is constant and satisfies the following conditions: 



r^o...^ ri23 77 = v , 



riTi V = T 2 T 2 r] = T 3 T 3 r] , 



rj = irj* . (A.52) 
Notice that, for the spinor e, the first of these projections can be recast in the form: 



r x o.. 



( cos aT 



f 123 ) 



123 - sinai 123 e = e . 



(A.53) 
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B Asymptotic behaviour of the fluctuations 



This appendix is devoted to the determination of the asymptotic form of the fluctuations 
around the kappa-symmetric static configurations of the D5-brane probe. In subsection B.l 
we will consider the large and small r behaviour of the solutions of the equation of motion 
corresponding to the lagrangian (7.5), which describes the small oscillations around the unit- 
winding embedding (7.2). In section B.2 we will study the asymptotic form of the fluctuations 
of the n-winding embedding in the abelian background. Following our general arguments, 
the UV behaviour of the abelian and non-abelian fluctuations must be the same and, thus, 
from this analysis we can get an idea of the nature of the small oscillations around the general 
non-abelian n-winding configurations (whose analytical form we have not determined) for 
large values of the radial coordinate. 

B.l Non abelian unit-winding embedding 

For large r the lagrangian (7.5) takes the form: 

C = - l -e^[r(d rX ) 2 + 2r X d rX + r X 2 + r (d xX ) 2 + r(d vX ) 2 } , (B.l) 

where we have not expanded the dilaton and we have eliminated the exponentially suppressed 
terms. Using the asymptotic value of d r (fi: 

d r <b w 1 — , (B.2) 

we obtain the following equation for £: 

^ + ^^+( M2 - i2 + 1 + f^:)« = - < B - 3 > 

To study the UV behaviour of the solutions of this differential equation it is interesting to 
rewrite it with the different coefficient functions expanded in powers of 1/r as: 



+ K + ^ + ^f + + ( & o + ^ + Jf - • • • K =0. (B.l) 



where 



a Q = 2 , b = M 2 - I 2 + 1 , ai = b 1 = \, a 2 = b 2 = ~ , ■ ■ ■ (B.5) 

Z o 

We want to solve eq. (B.4) by means of an asymptotic Frobenius expansion of the type 
£ = r p (c + Ci/r + • ■ •) for some exponent p. By substituting this expansion on the right- 
hand side of eq. (B.4) and comparing the terms with the different powers of r, we notice 
that there is only one term with r p , namely 6o c o rP > which cannot be canceled. In order to 
get rid of this term, let us define a new function w as 

t = e ar w, (B.6) 
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with a being a number to be determined. The equation satisfied by w is the same as that 
of £ with the changes: 



ai -> a, , bi-^bi + acii, (i = l,2, • • •) . (B.7) 

It is clear that we must impose the condition: 

a 2 + «a + 6 = , (B.8) 

which determines the values of a. Writing now w = r p (co + ci/r + ■ ■ •) and looking at the 

highest power of r in the equation of w (i.e. p — 1), we immediately obtain the value of p, 
namely: 

p = -^l±h, (B.9) 
and, clearly, as r — > 00 the asymptotic behaviour of £ (r) is: 

£( r ) ~ e Qr r"(l + o(i)). (B.10) 
In our case, it is easy to verify that the values of a and p are: 



a = -l±iVM 2 -P , / ? = ~4- ( B - U ) 
Then, it is clear that we have two independent behaviours for the real function £(r), namely: 

f(r) -cos(VM 2 -/ 2 r) , — sin ( VM 2 - Z 2 r ) . (B.12) 



Notice that all solutions decrease exponentially when r — > 00. 

Let us now turn to the analysis of the fluctuations for small values of the radial coordinate. 
Recall that r* < r < 00. Near r*, one can expand sin^o and cos^o as follows: 



sin^o ~ 1 — cothr* (r — r*) + (r — r*) 2 + 



11 + cosh r* 
2 sinh 2 r* 



y " coth. t \ 

cos# ~ \/2cothr*(r - r*) 1 — — (H -n — ) (r - r*) + 

v L 2 2 cosh r* 



. (B.13) 



Using these expressions it is straightforward to show that the lagrangian density of the 
quadratic fluctuations is given by: 

C = — e 3 *- [A(r)(d rX ) 2 + B(r)2xd r x + C(r)x 2 + D(r)(d xX ) 2 + E(r)(d v x) 2 }, (B.14) 
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where 0* = 4>( r *) and the functions A, B, C, D and E are of the form: 



A(r 

B(r 

C(r 

D(r 
E(r 



[8tanhr, (r -n) 3 ]2 (1 + A(r)) , 



2(r — r*) 
a/ coth r* 
1 



1 + B(r) ) , 



^2 cothr* (r - r*) 



1 +C{r)) , 



Vcothr* 



2(r-r*) (1 + V{r)) , 



[ tanhr* ^2(r - r») ( 1 + £(r) ) 



The functions „4, £>, C, X>, £, satisfy: 

A(r),B(r),C(r),V(r)£(r) ~ o(r - r*) . 



(B.15) 



(B.16) 



Notice that, remarkably, after integrating by parts the x^rX term in the lagrangian, the 
singular term of C(r) cancels against the leading term of B(r). The equation of motion for 
£ near r* is: 



p | A'(r) ^ | ^(r)-C(r)+M 2 D(r)-/ 2 E(r) ^ = Q 



A(r) 



A(r) 



(B.17) 



In order to solve this equation in a power series expansion in r - r„, it is important to 
understand the singularities of the different coefficients near r*. It is immediate that: 



A'(r) 



3 1 A'(r) 
+ 



+ regular . 



A(r) 2 r - r* 1 + A(r) 2 r - r* 
Similarly, the coefficient of £ has a simple pole near tv 

B'(r)-C(r) + M 2 D(r)-l 2 E(r) 3B '(r*) - C'(n) + 2M 2 r* - 2/ 2 tanhr, 1 



A(r) 



r — r* 



(B.18) 



+ regular 



(B. 



It follows that the point r = r* is a singular regular point. The corresponding Frobenius 
expansion reads: 



£(r) = (r-r») A ]T c n (r - r,)» , 



(B.20) 



n=0 



where A satisfies the indicial equation, which can be obtained by plugging the expansion in 
the equation and looking at the term with lowest power of r — r* (i.e. A — 2). In our case, A 
must be a root of the quadratic equation: 



A(A-1) + -A = , 



(B.21) 
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i.e.: 



A = 0,~. 



(B.22) 



This means that there are two independent solutions of the differential equation which can 
be represented by a Frobenius series around r = r*, one of then is regular as r — > r* (the one 
corresponding to A = 0), whereas the other diverges as (r — r*)~2. 



B.2 n- Winding embeddings in the abelian background 

Let us consider the abelian background and an embedding with winding number n, for which 
dpip = n, sin 9dg9 = nsin9 and — ip — constant. Let us define 



n (1 + cos#) n - (1 -cos#) n 1 

+ — cot u , 



n cos 6* 1 

= 2smT + 2 C ° 2sin# (l + cos#) n + (l-cos#) n ' 2 



nsin# 
W = - 



nan i 



v(n-l) 



2sin6> (1 + cos6>) n + (1 - cos6>) n 



(B.23) 



Choosing r as worldvolume coordinate and considering embeddings in which 9 depends on 
both r and on the unwrapped coordinates x, we obtain the following lagrangian: 



£ = - e 2 * sin 9 



( e 2h + y 2 + W 2j (J + ( e 2h + ^ 2 )((d r 0) 2 + (9 X 0) 2 )) + 



+ (e 2/l + W 2 )<9 r # - V 

We shall expand this lagrangian around a configuration #o( r ) sucn that 

2( r-r.) _ I (l + cosg (r))" + (1 - cos G (r)) n 



(B.24) 



(sin # (r)) 



n+l 



(B.25) 



Notice that #o( r ) remains unchanged under the transformation n — > — n. Thus, without loss 
of generality we shall restrict ourselves from now on to the case n > 0. Let us now define 



V (r) = V\ e{r )=e {r) , 
Using 



W (r) = W\ 



d r 9 



(r)=6 (r) , 
1 

" ~v ' 



<9V 



«(r)=fo(r) 



;B.26) 



(B.27) 



we obtain the following quadratic lagrangian 



C = 



D 2cf> 



sin^ (e 2ft + W$) 



e 2h + V 2 + W 2 

Vr 



(V^(d r x) 2 - 2V V oeX drX + 



+ fx 2 )+v (d xX f 



(B.28) 
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Keeping the leading terms for large r, the lagrangian becomes: 

C = -\<» [ ^±1 {drX f + 2 X d rX + x 2 + ^ fe) 2 ] , (B.29) 

and, if we represent x as i n ec L- (7-6) with / = 0, the equation of motion for £ becomes: 

1 / An 1 In 

^+(2 + -R£+ (M 2 + -—- 2 +—-- U = 0. (B.30) 

Now we have the following coefficients in eq. (B.4): 

An 1 1 

a = 2 , 6 = M 2 + , ai = - , &! = — — . (B.31) 

(n + iy 2 n + 1 

By plugging these values in eq. (B.8) we obtain the following result for the coefficient a of 
the exponential: 

'(^-j) 2 " ^ ■ (B.32) 



a = -1 ± 



that M 2 > (^i) 2 and define 



+ 

Let us distinguish two cases, depending on the sign inside the square root. Suppose first 

n-l ' " 
n+1' 

n — 1\2 
+ 

In this case the values of the exponent p obtained from (B.9) are: 



M 2 = M 2 - (^y) • (B.33) 



p = — - =F ^ 1 - % . (B.34) 
^ 4 ^ 2(n + l)M V ; 

It follows that the two real asymptotic solutions are: 

6 r ~ Tl — 1 6 r ~ Ti — 1 

f (r) cos [Mr - — — — — ^ logr] , — sin [Mr - logr] . (B.35) 

r* 2(n + l)M r* 2(n+l)M 

Both solutions vanish exponentially when r — > oo. 

If M 2 < 0, let us define M 2 = — M 2 . In this case a is real, namely a = — 1 ± M. Notice 
that M < 1 and thus a < 0. The independent asymptotic solutions are: 

f(r) ~ e (Ai ~ 1)r r^ 1 " T^TTTm) ? e -(*+i)r r -j (i+ (^tjm) ? (B.36) 

and both decrease exponentially, without oscillations, when r — > oo. This non-oscillatory 
character of the functions in eq. (B.35) make them inadequate for the type of boundary 
conditions we are imposing and, therefore, we shall discard them. 

Notice that, for n — 1, the large r asymptotic solutions (B.12) and (B.35) coincide. 
This is of course to be expected since the abelian and non-abelian configurations coincide in 
the UV. It is also interesting to compare the magnitude of the fluctuation with that of the 
unperturbed configuration for large r. By inspecting eq. (B.25) one readily concludes that 

9 (r) ~ e~^ r , (r -> oo) . (B.37) 
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By comparing this behaviour with eq. (B.35) one finds 



e "+ 1 



i 

f4 



(r — > oo) 



(B.38) 



0o(r) 

Thus, for n > 1 one has that — > as r — > oo. On the contrary for n = 0, both in the 

abelian and non-abelian case, the ratio -^Xj diverges in the UV and the first order expansion 
breaks down. 

Let us now consider the IR behaviour of the fluctuations. Near r* one has to leading 
order that sin 9 ps 1 and 



cos# 



VTT 



77/ 



n 



Vq ~ Vn 2 + 1 \Jr — r* , 

n 2 + 1 
K)0 ~ n — • 



(B.39) 



The IR lagrangian is of the same form as in eq. (B.14) (with E — since we are now 
considering the case in which x is independent of ip). The functions A(r) to D(r) are now 
of the form 



A(r) 
B(r) 
C{r) 



+ l)(r-r*) 2 ( 1 + o(r - r*)) , 



(n 2 + 1) 



V r - r * (l + o(r -r*)) , 



1 (n 2 + 1)2 



4 ^/r — r 
n 2 - ] 



( 1 + o(r - r* 



-D(r) = (r* H — ) V^ 2 + 1 \Jr - r\ ( 1 + o(r - r*)) . 



;b.4o) 



Notice that, also in this case, the coefficients of the functions above are such that, after a 
partial integration, the singular term of C(r) cancels against the leading term of B(r). The 
differential equation that follows for £ in the IR is: 



3 

2 r~ 



— + o(r-r*) W + o(^^)£ = , 

- r* ' r - r* 



(B.41) 



and, therefore, the indicial equation is the same as for eq. (B.17) (i.e. eq. (B.21)). It follows 
that also in this case there exists an independent solution which does not diverge when 
r — > r*. 
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